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1 .  Summary  of  Research  Objectives  and  Results 

In  this  research  program,  we  have  been  concerned  with  the  reactions  of  nitrogen  atoms  under 
upper-atmosphere  conditions.  Reactions  of  N(4S)  and  N*(2P,  2D)  with  molecular  oxygen  are  the 
major  source  of  vibrationally  excited  NO  in  the  upper  atmosphere,  which  in  turn  is  responsible  for 
infrared  emission  in  the  5  to  6  |im  region.  Measurement  of  the  NO  production  rates  and 
vibrational  distributions  are  carried  out  in  the  FACELIF  reactor  at  the  Air  Force  Geophysics 
Laboratory  (Hanscom  AFB,  Massachusetts)  using  resonance-enhanced  multiphoton  ionization 
(REMPI)  detection.  The  REMPI  detection  efficiency  for  N  atoms  and  NO  molecules  is  therefore 
of  considerable  interest. 

In  previous  work  on  this  system,  the  rates  and  branching  ratios  in  the  reaction 
N(4S)  +  02  ->  NO(v  <  7)  +  O 

w  ere  determined;  however,  an  unexpected  apparent  alternation  of  rates  into  neighboring  v  states 
required  further  investigation.  Similar  measurements  for  the  reaction 
N*  (2p,  2D)  +  02  -4  NO  (v  <  26)  +  O 

would  be  highly  desirable,  in  order  to  check  infrared  chemiluminescence  results  from  the 
COCHISE  reactor.  In  order  to  optimize  such  an  experiment,  rates  for  deactivation  of  the 
metastable  nitrogen  atoms  by  constituents  in  the  flow  reactor  must  be  known.  During  the  course  of 
this  research,  a  report  appeared  from  the  SRI  group  [J.  Phys.  Chem.  22,  5977  (1988)]  in  which 
the  rate  coefficient  for  the  N(2D)  +  0(3P)  quenching  reaction  was  reported  as  ca.2  x  1(H 1  cm3sec'1 
at  300K,  considerably  larger  than  previous  measurements.  We  undertook  to  remeasure  this 
reaction  rate  in  the  FACELIF  reactor.  The  0(3P)  atoms  were  generated  by  titrating  NO  directly 
into  the  active  nitrogen  flow.  Analysis  of  the  results  required  that  the  rate  coefficient  for  the 
reaction  N(2D)  +  NO  — »  N2  +  O  be  known  accurately,  and  this  was  also  determined.  A  finite 
mix'ng  time  correction  is  also  necessary.  The  best  estimate  of  the  rate  coefficient  from  this  work  is 
(6.9  +  7,  -1.1)  x  10 '13  cm3sec'1  at  T  =  298  K,  considerably  smaller  than  the  rate  reported  by  the 
SRI  group  and  in  much  better  accord  with  values  required  by  atmospheric  models.  These  results, 
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described  in  Section  2  of  this  report,  are  in  the  process  of  publication  in  the  Journal  of  Chemical 
Physics,  and  complete  kinetic  data  are  documented  in  AFGL  technical  report. 

The  mixing-time  and  other  corrections  indicated  by  the  work  on  the  N(2D)  +  0  quenching 
rate  suggested  that  a  previous  report  by  Phillips  et  al.  [J.  Phys.  Chem.  21, 5001  (1987)],  in  which 
the  N(2P)  +  O2  reaction  rate  was  determined  to  be  (1.8  ±  0.2)  x  10*12  cm3  sec*1,  might  also  be  in 
error.  Preliminary  results  of  a  re-measurement  of  this  rate  in  the  FACELIF  reactor,  described  in 
Section  3,  give  a  value  of  ~3.4  x  10*12  cm3sec*1,  in  much  better  accord  with  earlier  measurements 
by  Kaufmann  and  co-workers  [J.  Chem.  Phys.  21,  4701  (1980)].  However,  an  independent 
measurement  by  Piper  et  al.  [private  communication)  gives  a  value  of  k[N(2P)  +  O2]  ~  2.  x  10*12 
cm3sec*1,  which  agrees  with  the  previous  Phillips  et  al.  measurement.  This  issue  is  unresolved  at 
the  present  time. 

During  the  course  of  this  work,  we  observed  that  N(4S)  atoms  could  be  efficiently  detected 
by  Resonantly  Enhanced  Multiphoton  Ionization  (REMPI).  The  magnitude  of  the  ion  signal,  after 
correction  for  laser  power  dependence,  is  directly  proportional  to  the  N(4S)  concentration  and, 
once  calibrated,  is  a  direct  and  highly  sensitive  method  for  N(4S)  concentration  measurements  in 
fast-flow  kinetics  experiments.  These  results,  described  in  Section  4,  are  in  the  process  of 
publication  as  a  Note  in  Spectrochimica  Acta. 

It  was  noted  earlier  that  the  vibrational  state  distribution  in  the  N(4S)  +  O2  reaction,  as  probed 
by  REMPI,  showed  an  unexpected  alternation  between  adjacent  vibrational  states.  It  was 
suggested  by  C.  Feigerle  [private  communication]  that  this  could  be  due  to  partial  saturation  of  the 
initial  (A  <-  X)  absorption  step  in  NO.  Accordingly,  a  detailed  kinetic  and  spectroscopic  model  for 
(1  +  1)  multiphoton  ionization  of  Hund's  case  (a)  diatomic  molecules  in  this  regime  was 
investigated.  The  model  predictions,  verified  for  NO(v=l)/NO(v=0)  signal  ratios  and  for  REMPI 

fWpr-tinn  nf  hitrVi  yihratirmal  nf  M->  CYlTVt  in  nitrnoen  Hicr,tiarar*Q  cnoapct  that  thp  nrf*v’Ol,<; 

measurements  of  NO  vibrational  state  distributions  from  the  N(4S)  +  O2  reaction  carried  out  in  the 
FACELIF  reactor  do  indeed  need  to  be  corrected  for  partial  saturation  effects  to  obtain  accurate 
values  of  k(v).  These  results,  described  in  Section  5,  are  being  prepared  for  eventual  publication. 
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In  order  to  completely  reduce  the  data  from  FACELIF  experiments  on  (N  +  O2)  or  (N*  + 

-» NO  (v),  it  also  necessary  to  know  vibrational  deactivation  efficiencies  for  highly  vibrationally 
excited  levels  of  NO.  We  have  investigated  the  feasibility  of  using  Stimulated-Emission  Pumping 
(SEP)  to  generate  the  high  instantaneous  densities  of  NO(v")  required  for  such  measurements. 
While  SEP  via  the  B(2n)  state  using  a  commercially  available  ArF  excimer  was  found  not  to  be 
feasible,  a  two-photon-pump  scheme  via  the  A2!  state,  followed  by  stimulated  A  ->  X  emission, 
did  succeed  in  moving  population  into  v"  =  4  levels  of  the  X  state,  as  described  in  Section  6. 
However,  due  to  limitations  on  the  availability  of  laser  sources,  it  was  not  possible  to  carry  out 
detailed  measurements  on  NO  (v")  state  relaxation.  In  a  concluding  section,  several  possible 
directions  for  future  investigations  are  indicated. 
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2.  Quenching  of  N(2D)  by  0(3P) 

[for  publication  in  the  Journal  of  Chemical  Physics,  v.  92  (April  15  1990)] 
1 .  Introduction 


N(  D)  is  an  important  constituent  of  the  mesophere  and 

1-3 


thermosphere.  Atmospheric  modeling  studies 
the  reaction 


indicate  that 


N{2D)  +  02  -  kfo->  NO { v)  +  0(3P,1D) 


(1) 


is  the  primary  source  of  the  NO  concentrations  first  observed  in 

4 

rocket  experiments  .  While  the  value  of  the  rate  coefficient 
for  reaction  (1)  is  well  established^-^,  the  coefficient  for 
the  cuenching  reaction 


N(2D)  +  0(3P)  N(4S)  +  0(3P) 


(2) 


is  still  the  subject  of  some  dispute.  The  relative  magnitude  of 
these  two  rate  coefficients  is  important  in  determining  the 
partitioning  between  the  two  product  channels.  Modeling 

calculations  require  that  reaction  (1)  be  the  dominant  loss 

2  11 
channel  for  N(  D) .  A  recent  experimental  measurement  of  the  rate 

-11 

coefficient  for  reaction  (2)  reported  a  value  of  3.4  x  10  x 

3  -1  -11 

exp(-145/T)  cm  sec  ,  which  would  correspond  to  2.1  x  10 

3  -1 

cm  sec  at  T  =  300  K;  the  quoted  error  limits  are  on  the  order  of 
±50%.  This  value  is  considerably  higher  than  previous  determinations 
and  more  discordant  with  preferred  modeling  values2' 
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•  -12  3  -1 

which  require  kQ  <  10  cm  sec  .  If  this  rate  were  correct,  then 
2 

N(  D)  would  contribute  to  the  overall  depletion  of  NO  by 

,  4 

enhancing  the  N(  S)  concentration  which  reacts  rapidly  with  NO 
through  the  reaction 

N(4S)  +  NO  N2  +  0  (3P)  (3) 

.  2 
This  would  necessitate  a  re-evaluation  of  the  role  of  N (  D)  in 

the  upper  atmosphere  and  the  identification  of  another  source 

term  for  NO. 

We  have  re-measured  the  rate  coefficients  for  both  reactions  (1) 
and  (2)  in  order  to  establish  the  true  partitioning  between  these  two 

2  i  r 

N (  D)  loss  channels.  This  was  done  using  the  FACELIF  reactor  with 

multi-photon  ionization  (MPI)  detection  of  reactant  and  product 

species.  The  experimental  method  makes  use  of  reaction 

3  11 

(3)  to  generate  the  0{  P)  atom  quenchers  .  As  this  reaction 
takes  place  directly  in  the  main  flow,  recovery  of  kQ  requires 

that  the  rate  coefficient  for  the  reaction 

N(2D)  +  NO  —  kNQ..  >  N2  +  o  (4) 


be  known  accurately.  Reported  values17'18  for  this  rate  coefficient 

range  from  5.9x10  ^  to  1.8xl0-10  cm3s-1  ;  thus,  it  was  necessary  to 

carry  out  a  measurement  of  k.„  as  well,  in  order  to  determine  the 

nO 

0‘ 


correct  value  of  k 
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2 .  Experimental 

2.1.  Apparatus 

The  experimental  set-up  was  essentially  the  same  as  that  used  in 

1 6 

previous  rate  measurements  performed  using  the  FACELIF  reactor  and 

will  therefore  be  described  only  briefly  here.  A  schematic  is  shown 
.  .  2 

m  Fig.l.  The  N(  D)  atoms  were  produced  in  a  70W  microwave  discharge 
of  a  1-2%  mixture  of  ^  in  helium.  The  N,,  and  He  were  purified  using 

gettering  furnaces  containing  Cu  at  6U0°  C  and  Ti  at  800°  C, 
respectively,  to  remove  00  impurities  prior  to  entering 

the  discharge  tube.  The  NO  (Matheson  98%  purity)  was  flowed 
through  a  purifying  manifold  consisting  of  an  ascarite  trap  and  a 
liquid  nitrogen/methanol  slurry,  and  was  then  diluted  with  He  to 
-10%  concentration  prior  to  use.  All  other  quenchers  were  used 
without  furthur  purification. 

Metastable  nitrogen  atoms  are  deactivated  at  essentially  every 

wall  collision^'  19 .  The  N<^D)  concentrations  created  by  microwave 

discharges  are  approximately  two  orders  of  magnitude  lower  than  the 

ground  state  nitrogen  atom  concentration.  The  combination  of  these 

two  factors  made  it  essential,  under  certain  circumstances,  to  be 

able  to  introduce  the  metastable  atoms  into  the  reaction  zone  as 

2 

quickly  as  possible,  while  the  N(  D)  was  still  in  sufficient 
concentration  to  provide  an  adequate  signal  level.  This  necessitated 
a  fast  pumping  speed  which  was  provided  by  a  Aoots  blower  (Leybold- 
Heraeus  WS1000)  backed  by  a  large  displacement  forepump  (Heraeus- 
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Englehaird  DK180)  .  The  maximum  line*,  i  ow  velocity  attainable  usir.j 

this  configuration  was  of  che  order  of  6000  cm  s-1,  which  correspond- 

to  a  minimum  flow  time  of  approximately  2  ms.  Longer  flow  times  could 

be  achieved  by  reducing  the  pumping  speed  as  required.  Quenching 

gases  were  injected  through  a  loop  injector  whose  position  could  be 

adjusted,  over  a  limited  range,  along  the  length  of  the  flow  tube, 

again  altering  the  flow  time.  Some  cnparative  studies  were 

performed  using  a  perforated  ball  injector  to  probe  any  possible 

disparities  in  mixing  properties.  These  indicated  that  the  mixing 

characteristics  were  comparable  with  marginally  poorer  mixing  for 

the  ball  injector. 

2 

The  N (  D)  concentration  was  monitored  using  resonance  enhanced 

18  19 

multiphoton  ionization  '  .  The  frequency  doubled  output  of  a 

Nd : YAG  pumped  dye  laser  (Quanta-Ray  DCR2A,  PDL-1,  WEX-1)  was  used  to 
2 

photo-ionize  N(  D),  resonantly  enhanced  by  the  two-photon 

.  .  20 
transitions  at  268.95  and  268.98  nm  .  The  current  produced 

across  a  pair  of  Ni  wire  grids,  biased  at  90V,  was  amplified 

using  an  Ithaco  model  1211  current  preamplifier  and  averaged  on 

an  EG  &  G  model  162  Boxcar  Averager.  Signal  levels  were  generally 

measured  by  setting  the  laser  on  resonance  at  each  data  point  and 

subtracting  the  off-resonance  baseline  from  the  averaged  signal. 

Some  results  were  also  taken  by  scanning  the  laser  across  the 

resonance  lines  with  no  significant  difference  in  the  results. 

3 

The  method  used  to  generate  the  Q(  P)  quenchers  involves  the 

titration  of  NO  directly  into  the  f lowtube'*'^ .  The  NO/He  mixture  is 

introduced  into  the  flow,  through  the  loop  injector,  where  it  reacts 
.  .  4  21 

rapidly  with  N{  S)  ,  which  is  present  in  high  concentrations  in 
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4 

microwc.ve  discharge  flows  of  ^/He.  Under  ideal  conditions,  the  N(  S) 

will  oe  in  sufficient  excess  that  all  the  NO  will  be  rapidly 
3 

converted  to  0(  P)  before  making  any  significant  contribution  to  the 
2 

N (  D)  quenching.  The  diffusion  coefficient  for  0  atoms  in  He  is  700 
cm2  sec-1  22,  and  the  surface  recombination  coefficient  is  on  the 
order  of  10-<1  22 .  For  the  flowspeeds  and  pressures  used  in  our  2  in. 
diameter  flowtube,  we  estimate  less  than  1%  loss  of  0(3P)  over  the 
reaction  length. 

2 . 2  Experimental  Procedure 

Preliminary  experiments  were  carried  out  investigating  the 

2 

variation  in  the  total  N(  D)  quenching  rate  over  a  wide  range  of 
reaction  times.  These  were  used  to  calculate  an  effective  rate 
coefficient  at  each  reaction  time.  This  is  defined  as 

keff  =  -  a(v/d) (dln[N(2D) ]/d[Q])  (5) 

where  v  is  the  average  flow  speed,  d  is  the  injector  to  detector 

distance,  a  is  the  non-laminar  flow  correction  having  a  numerical 
q 

value  of  1.62  ,  and  [Q]  is  the  initial  quencher  concentration.  The 

use  of  Eq.  (5)  is  equivalent  to  assuming  that  all  the  initially  added 

3  4 

NO  is  converted  instantaneously  to  0(  P)  by  the  N(  S)  without 

2 

contributing  to  the  N(  D)  decay.  The  results,  shown  in  Fig. 2, 
display  a  rapid  decrease  in  the  effective  rate  coefficient  with 
increasing  reaction  time.  Even  at  ~50ms  there  is  some  indication 
that  the  rate  coefficient  is  still  approaching  an  asymptotic 
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-1°  3  -1 

value  of  less  than  1x10  “  cm  s  .  l’he  flow  conditions  over  this 

range  of  reaction  times  necessarily  varied  by  a  substantial  amount. 
The  total  pressure  varied  from  0.2  to  4  Torr  in  going  from  short  to 

4 

longer  reaction  times.  There  was  a  corresponding  change  in  the  N(  S) 

concentration,  although  this  was  not  explicitly  measured  for  all 

points  along  the  curve.  The  percentage  of  initially  added  NO 

remaining  unreacted  at  the  detection  point  was  also  measured  for  the 

24  25 

same  range  of  reaction  times.  The  NO  y(0,0)  MPI  signals  '  were 
used  for  this  purpose.  The  magnitude  of  this  signal  is  directly 
proportional  to  the  NO  concentration.  A  blank  titration  was  performed 
under  the  same  set  of  experimental  conditions,  except  with  the 
discharge  turned  off.  The  ratio  of  the  two  slopes  yields  a  value  of 
the  fraction  of  initially  added  NO  that  remains  unreacted.  The 
results  are  shown  in  Fig. 3. 

Comparison  of  Figs.  2  and  3  shows  clearly  that  the  faster 
effective  rate  constants  are  correlated  with  significant  residual 

NO  concentrations.  The  short  reaction  time  measurements  represent 

.  .  .  2 
a  region  in  which  reaction  with  NO  dominates  the  N(  D)  decay. 

These  measurements  indicate  an  extreme  lower  limit  on  k.,^.  The 

NO 

longer  reaction  time  measurements,  to  the  extent  that  there  is 

still  some  NO  contribution  to  the  effective  decay  rate,  represent 

3 

an  extreme  upper  limit  on  the  0(  P)  quenching  rate  coefficient. 

The  limiting  values  obtained  for  kNQ  and  kQ  obtained  in  this  way, 

indicate  a  ratio  for  these  two  rate  coefficients  in  excess  of  25 
to  1.  Consequently,  for  short  reaction  times,  where  [NO]  ~  [0],  the 

3 

0(  P)  contribution  is  smaller  by  a  factor  k^/k^,  and  may  be 
legitimately  neglected  when  determining  k^.  Conversely,  the 
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large  value  of  kNQ  compared  to  kQ  means  that  even  very  small 

amounts  of  NO  may  cause  a  significant  degree  of  quenching. 

It  has  been  reported  that  REMPI  detection  of  N(2D)  at  269  nm  can 
2 

generate  N(  u) through  the  two  photon  dissociation  of  NO  26,27.  jn 

order  to  ensure  that  this  effect  was  not  compromising  our  results  we 

o 

looked  for  the  characteristic  N(  D)  signal  at  high  NO  flow  with  the 
discharge  turned  off  (i.e.  with  [NO]  significantly  higher  than  normal 
operating  conditions) .  No  signal  was  observed  until  [NO]  was  over 
100  times  that  used  during  titrations.  The  laser  power,  measured 
after  exiting  the  flow  tube,  was  <0.5  mJ  pulse-1  compared  to  the  2-5 
mJ  pulse-1  reported  by  Jusinski  et  al.  26 .  This  gives  a  production 
rate  ratio  of  between  1:16  and  1:100  for  the  two  systems.  Also,  NO 
concentrations  in  the  detection  region  were  always  <0.5  mTorr.  The 
combination  of  these  two  factors  accounts  for  the  lack  of  any 
measurable  levels  of  N(2D)  production  from  NO  dissociation. 

4 

Measurements  of  kNQ  and  kQ  both  require  knowledge  of  [N (  S) ]  to 

determine  how  fast  the  NO  is  depleted  and  therefore  the  extent  to 
.  .  2 

which  it  contributes  to  the  N(  D)  decay.  For  long  reaction  times 

28  4 

the  standard  NO  titration  method  may  be  used  to  determine  [N (  S) ] 
accurately.  This  technique  has  serious  limitations  at  the 
shorter  reaction  times  to  which  the  NO  rate  coefficient 
measurements  were  constrained.  To  obtain  a  distinct  titration  end 
point,  reaction  (3)  must  have  sufficient  time  to  go  to 
completion.  The  method  therefore,  is  best  suited  to  conditions 

4 

where  [N {  S) ]  is  high  and  the  reaction  time  is  long.  This  is 
precisely  the  opposite  of  the  requirements  for  accurate  kN(^ 
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measurements.  Other  practical  difficulties  were  encountered  in 

the  form  of  greatly  reduced  signal  levels  and  interfering 

fluorescence.  The  fast  flows  and  consequent  low  pressures  (0.2- 

0.4  Torr)  lead  to  significantly  reduced  signal  levels  and  a 

subsequent  reduction  in  signal  to  noise.  In  addition,  stray  light 

from  the  discharge  created  a  background  large  in  comparison  to 

the  target  signal.  Anomalous  effects  were  observed  for  small  NO 

additions  in  the  form  of  an  initial  rise  in  signal  level.  It  is 

possible  that  this  has  its  origins  in  fluorescence  resulting  from 

29 

some  coincident  reaction.  It  has  been  previously  reported 

4 

that  some  short-lived  source  of  N£ (B) ,  other  than  N(  S)  atom 

recombination,  exists  in  discharge  flows  which  may  interfere  with 
titration  measurements. 

Attempts  were  made  to  fit  the  titration  plots  for  short  reaction 

times.  This  introduced  a  further  limitation  beyond  the  practical 

problems  already  mentioned.  Fitting  would  require  the  use  of  the  rate 

coefficient  for  reaction  (3),  which  has  associated  error  bounds  that 
.  4 

will  propagate  into  the  (N(  S) )  value.  As  a  result  of  these 
limitations,  an  alternative  method  of  determining  [N (  3)]  was 
devised.  Measurements  of  NO  depletion  were  made  using 

J  (0,0)  MPI  band  intensities,  as  described  earlier.  The  high 
sensitivity  of  MPI  to  NO  enables  very  small  concentrations  to  be 
detected  and  it  is  therefore  possible  to  work  with  [N (4S) ) » (NO) .  The 
NO  decay  can  then  be  described  by  the  p:. iudo  first-order  rate 
equation 


[NO]  =  [NO] Q  exp(-kN[N(4S)  ]t) 


(6) 
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4 

where  [N0]q  is  the  initial  concentration  of  NO  added.  The  N(  S) 
concentration  is  simply  given  by 

[N(4S)]  =  -ln{[NO]/[NO]Q}/(kNt)  (7) 

Since  the  ratio  [NO]/[NO]o  may  be  determined  very  precisely,  the 
accuracy  of  this  technique  is  limited  by  the  uncertainty  in  k^. 

For  slow  flows  (  >  20ms),  it  is  reasonable  to  assume  that  the 
mixing  time  is  short  in  comparison  to  the  total  flowtime.  For 
very  fast  flows  (  <  5ms) ,  this  assumption  may  lead  to  a  significant 
underestimation  of  the  rate  coefficients.  A  standard  technique 
used  to  probe  mixing  effects  is  to  perform  rate  measurements  over 
a  range  of  injector  positions  under  constant  flow  conditions.  A 
plot  of  the  decay  coefficient  as  a  function  of  the  flowtime  will 
give  a  straight  line  whose  x-axis  intercept  gives  a  value  for  the 
effective  mixing  time  and  whose  slope  yields  a  value  for  the  rate 
coefficient.  In  the  case  of  NO,  the  accompanying  reaction  with 

4 

N (  S)  depletes  NO  and  thus  compromises  the  technique.  In  order  to 

circumvent  this  problem,  we  performed  the  rate  measurements 
substituting  O2  as  a  quencher.  If  it  is  assumed  that  the  mixing 

characteristics  for  O2  and  NO  are  comparable,  then  the  effective 

mixing  time  obtained  using  O2  may  be  applied  to  NO,  provided  that 

identical  flow  conditions  are  employed. 
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3.  Analysis  of  Results 

Measurements  of  both  and  k^  involve  injecting  NO  molecules 

directly  into  the  main  discharge  flow,  where  they  react  rapidly 
•  4  3 

with  N(  S)  atoms  to  generate  0(  P) .  Whatever  the  flow  conditions, 

the  NO  will  be  present  for  some  finite  time  in  the  tube.  The 
2  3 

N(  D)  atoms  may  be  quenched  by  0(  P)  and  NO,  both  of  which  will  be 

present  in  varying  concentrations  along  the  length  of  the 

flowtube.  From  reactions  (2)  and  (3),  the  rate  equation 

2 

describing  the  decay  of  N!  D)  is 

-d[N(2D)]/dt  =  kNQ[N(2D) ] [NO]  +  kQ[N(2D)][0]  +  kw[N(2D)]  (8) 

where  kw  is  the  first  order  rate  coefficient  for  wall  deactivation. 
Integration  of  this  equation  over  the  reaction  time  gives 

-ln{  [N(2D)  ]  /  [N  <  2D)  3  0 }  =  kNQ  Jo  [NO]  dt  +  kQ  fj  [0]dt  +  kwt  (9) 

where  [N(2D)]q  is  the  initial  concentration  of  N(2D).  The 
.  3 

concentration  of  0(  P)  at  any  given  time  is  simply  equal  to  [NO]o  - 
[NO] t,  the  amount  of  NO  reacted,  if  there  is  no  other  significant 
loss  of  0(3P)  . 

The  rate  of  change  of  NO  with  time  is  given  by 
-d [NO] / dt  =  kN[NO] (N(4S) ]  +  kNQ[NO][N*] 


(10) 
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where  [N  ]  is  the  metastable  nitrogen  atom  concentration.  Since 
4  * 

[N (  S)]»[N  ],  the  second  term  in  equation  (10)  may  be  neglected, 
and  the  NO  concentration  is  given  by  Eq.  (6) .  Substituting  this  into 
Eq.  (9)  gives 

-ln{  [N(2D)  ]/[N'(2D)  ]q}  =  kQ  [NO]  0t  +  (k^Q  -  kQ)  /kN  [N  (4S)  ] 

x  {1  -  exp (-kN [N (4S) ] t) } [NO] Q  (11) 


.  ,  4 

For  the  situation  where  [N(  S) ]  cannot  be  considered  to  be  in 
large  excess,  a  slightly  more  complex  expression  results.  The  NO 
concentration,  in  this  instance,  is  expressed  as 


-d[NO]/dt  =  V  [NO]o-[NO]t  }{[N(4s))0  -  (12) 


4 

where  [N(  S)]Q  is  the  initial  ground  state  nitrogen  atom 

concentration.  This  is  a  second  order  rate  equation  which  when 
integrated  gives*’ 


([NO]„-[n(4s)]) 


N4s)]  ([NO]0  —  [NO], ) 


[N°]o  |  n(4s)  -[NO],) 


=  ktft 


(13) 


Re-arranging  this  expression  for  [NO]  gives 
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[NO],- 


[N(4s)].[NO]oIgle>>0-H4#N' 


[n(4s)]  -([NO]0-[N(4s)])^,r 
[NO]0  * 


(14) 


2 

The  N (  D)  concentration  at  any  point  along  the  flowtube  is  given 
by 


=  km[NO}0t + (k0  -  *No)Jo[N°]dt + const- 


(15) 


J0maX  (NO]dt  must  be  evaluated  by  numerical  integration  over  the 

range  t  =  0  to  t  where  t  is  the  reaction  time  corrected  for 
max  max 


mixing  effects. 


4.  Results  and  Discussion 


An  estimate  of  the  effective  mixing  time  for  the  NO  kinetic  runs 
was  obtained  from  the  0 ^  rate  measurements.  The  results  are  shown  in 

Fig.  4.  The  total  pressure  in  the  flowtube  was  ~0.4  Torr  with  a 
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linear  flow  velocity  of  ~6000  cm  s  ^ .  A  weighted  linear  regression 
fit  to  the  data  yields  a  value  of  (0.48  ±  0.19) ms  for  the  effective 

mixing  time.  The  slope  of  the  line  in  Fig. 4  gives  a  value  of  the 
N(^D)  +  0 2  rate  coefficient  of  (5.89  ±  0.40)  x 

_ Y2  2  —1  5““10 

10  cms  ,  in  good  agreement  with  previous  determinations 

The  NO  experiments  were  carried  out  using  identical  flow 

4 

conditions  t.o  those  for  the  0^  runs.  The  N(  S)  concentrations 

were  varied  between  0.25  and  0.45  mTorr.  For  these  conditions 
Jo  (OJdt  -  Jo  [NO]dt,  but  as  kNQ  »  ,  the  0(~*P)  atom  quenching 

2 

may  be  neglected.  From  Eq.  (9),  a  plot  of  the  N(  D)  signal  decay 
versus  Jo  [NO]dt  will  give  a  straight  line  whose  slope  gives  a  value 
for  kNQ  directly.  Such  a  plot  is  shown  in  Fig. 5  for  t  =  3.90  ms  and 

[N(4S) ]  =0.315  mTorr. 

There  are  a  number  of  sources  of  error  that  must  be  considered 
in  calculating  kN0.  Uncertainties  in  signal  levels  are  shown  as 

.  2 
error  bars  in  Fig.  5.  The  dynamic  range  of  the  N(  D)  decay  which 

could  be  followed  was  limited  by  the  appearance  of  the  NO  7(0,4) 

2 

band  at  higher  NO  concentrations,  which  obscured  the  N(  D)  signal 
to  some  extent.  This  baseline  variation  is  reflected  in  the 
larger  error  bars  for  high  [NO]  additions.  The  vibrationally 

30 

excited  NO  population  is  presumably  the  result  of  the  reaction 

N2(A)  +  NO (X)  — >  N2(X)  +  NO (A)  (16) 

followed  by  y-band  emission,  although  this  has  not  beer,  confirmed 
experimentally . 
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The  vaL'ue  used  for  k„  in  the  calculations  is  taken  from  the 

N 

21 

literature  an <d  has  an  associated  uncex.  -'f  +26%.  This 

value  is  use  <3  j.ri  calculating  [N(4S)]  and  aga— .  the  integration 

of  lo  [ll0)cit,  This  error  manifests  itself  as  a  shift  in  the  x-axis 
position  of  "the  data  points.  The  extreme  values  for  k^  ,  [N(^S))f 

and  wer-eused  to  calculate  upper  and  lower  bounds  for  k^Q. 

These  bobnds  are  shown  in  Fig.  6  for  the  data  shown  in  Fig.  5. 

The  upp^a:  and  lower  bounds,  together  with  rate  coefficient 

obtained  usi-tg  the  mean  values  of  k.,  ,  (N(^S)l  and  t  ,  are  given 

N  max 

in  Tabl^  1  ffot  all  kinetic  runs.  The  error  limits  on  each  value 
are  the  statistical,  errors  {  +1  O  )  obtained  from  a  weighted  linear 
regression  £: it  "to  the  data  points  in  each  limit.  These  are  small 

in  compgUriscon  t  o  the  error  bounds  indicating  that  the  dominant 
source  o f  errror  arises  from  the  uncertainty  in  the  values  of  k^ 

and  the  hiixx-ng  time  .  The  final  values  are  obtained  by  taking  the 

weighted  mean  of  all  kinetic  runs.  This  procedure  yields  a  result 

of  (6.70,  +X.36,  -1.21)  x  10_11cm3s''1  for  kM^. 

NO 

The  ideal  conditions  in  which  to  measure  k^  are  at  long  reaction 

times  with  C^o]  ^(M  (^S)  ]  .  If  these  conditions  can  be  fulfilled  then 
Eq.  (11)  irax/  be  used  to  evaluate  k^.  Given  the  great 

disparity  between  kNQ  and  k^,  these  conditions  are  extremely 

difficuit  to  realize  in  practice.  At  the  slowest  flows,  the  use 

of  high  flov?tut>®  pressures  was  not  sufficient  to  offset  the  loss 

of  signal  dtJC  to  increased  wall  deactivation,  thus  constraining 

the  experiments  1  range  of  t.  The  maximum  [N(^S>]  values  obtained, 

even  under  tdheSC  extreme  conditions,  were  1.6x10^  cm  3  (~  5mTorr)  . 
The  low  valmc  c>i  k^  required  NO  additions  of  up  to  3.5mTorr  to 
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give  [N (2D) ] / [N <2D) J  <  0.05  and  thus  the  [N (4S) ] » [NO]  condition 

could  not  be  met.  Some  attempts  were  made  to  enhance  N-atom 

31 

dissociation  by  adding  SFg  to  the  discharge  flow  .This  did 

4 

result  in  a  significant  increase  in  [N  (  S) ]  (by  a  factor  of  up  to 

? 

2.5)  and  also  m  the  N(‘D)  signal  level.  Unfortunately,  the  use 
of  SFg  also  resulted  in  the  appearance  of  a  large  molecular  MPI 

spectrum  in  the  same  frequency  range  which  also  varied  with 

quencher  concentration  and  therefore  compromised  the  results.  It 

is  possible  that  such  interferences  would  not  occur  for  other 
2 

N (  D)  detection  methods.  If  this  were  true,  then  this  method  could 

4  2 

be  used  to  enhance  both  N(  S)  and  N(  D),  assuming  that  the 
species  generated  by  discharging  SFg  do  not  participate  in  the 

chemistry  in  any  way. 

Eq. (15)  was  used  in  calculating  kQ.  The  previously 

determined  value  of  k.,„  was  used  and  the  NO  contribution  to  the 

NO 

decay  subtracted  at  each  point.  Figs.  7  and  8  show  the  effect  of 
applying  this  correction  to  the  raw  data  for  two  separate  kinetic 
runs.  The  relative  contribution  from  NO  increases  with  increased 
quencher  addition  and  shorter  reaction  time.  Even  for  the  highest 

4 

(N (  S) ]  values  and  longest  flowtimes,  NO  accounted  for  ~30%  of 

9 

the  N(“D)  decay. 

4 

The  N(  S)  concentration  can  be  measured  by  the  more  usual 

4 

titration  technique  for  these  slow  flow  conditions.  [N(  S) ]  was 
measured  before  and  after  each  kinetic  run  with  results  agreeing 
to  within  2%  in  all  instances.  We  conservatively  estimate  an 
accuracy  of  +5%  for  the  NO  titration  technique.  The  mixing  time 
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was  not  explicitly  measured  for  these  longer  reaction  times  as 
this  can  reasonably  be  expected  to  be  short  compared  to  the 
flowtime.  Given  the  long  flowtime  and  large  N(^S)  concentration, 
lo  [0]dt  -4  [NO]Qt  and  the  correction  for  the  NO  contribution  to 

the  decay  involves  primarily  a  y-axis  adjustment.  The  magnitude 
of  this  adjustment  is  sensitive  to  the  exact  values  used  for 

and  k^.  *n  the  interests  of  self-consistency  within  the 

analysis,  when  one  extreme  value  of  kv,  was  taken  in  the 

N 

error  calculation,  the  value  of  k^  calculated  using  that  same 

value  was  used.  Thus  the  k„  and  k,,,.  errors  are  correlated  and,  in 

this  instance,  tend  to  cancel  one  another.  The  measurement  of  is 

therefore  less  sensitive  to  the  uncertainty  in  kQ  than  is  kNQ.  The 

error  bounds  for  k^  were  calculated  using  the  extreme  values  of 
4 

[N  (  S)  ] ,  kN,  and  kIJQ.  The  results  are  shown  in  Figs.  9  and  10  for 

the  data  in  Figs.  7  and  8,  respectively.  The  results  of  twelve 
separate  kinetic  runs  are  given  in  Table  2.  As  before,  the  errors 
quoted  for  each  value  are  the  statistical  values  in  each  limit 
and  are  small  compared  to  the  uncertainities  arising  from  the 

correlated  systematic  errors.  The  weighted  averages  of  these  data 
give  a  value  for  kQ  of  (6.93  +  0.71,  -1.08)xl0-13  cm3s-1.  The 

absence  of  any  systematic  tariation  in  kc  with  t  confirms  that  there 
is  no  significant  loss  of  along  the  flow  tube.  The  O2  formed 

would  react  rapidly  »-ith  \  (2D),  giving  an  apparent  increase  in  ko  at 
longer  reaction  times. 
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5.  Conclusions 

The  previously  reported  experimental  and  modeling  values  of  kQ 
are  listed  in  Table  3.  The  value  of  kQ  reported  here  is  in  good 

agreement  with  the  most  recent  experimental  measurement  made  by 

32  3  2 

Piper  ,  who  utilized  0(  P)  atom  generation  and  N(  D)  detection 

techniques  differing  from  those  employed  here.  It  is  also  in  fair 

33 

agreement  with  the  earlier  results  of  Davenport  et  al.  and 

9 

Ianuzzi  et  al.  .  The  major  conflict  is  with  the  result  of 
11 

Jusmski  et  al.  ,  who  recently  rekindled  the  ongoing  controversy 

regarding  the  disparity  between  experimental  measurements  and  the 

preferred  modeling  values  by  suggesting  that  the  true  rate  was  of 

-11  3  -1 

the  order  of  2.1x10  cm  s  .  The  experiment  by  Jusinski  et  al .  and 

the  studies  reported  here  were  performed  using  an  essentially 

identical  experimental  set-up  and  methodology.  The  origin  of  the 

disparity  in  the  results  is  not  immediately  apparent  although  it  is 

interesting  to  note  that  our  preliminary  investigations  (as 
illustrated  in  Fig.  2)  also  yield  a  value  of  kgf^  ~  2xl0~H 

3  -1 

cm  sec  at  t  =  7.5ms.  Our  results  and  calculations  indicate  that, 

4 

for  the  conditions  reported  in  ref.  11  ([N(  S) ]  =  lOmTorr,  t  =  7.5 

ms) ,  an  effective  rate  constant  (which  is  what  was  actually  reported) 
-12  3  -1 

of  -  2x10  cm  s  should  be  observed.  The  order  of  magnitude 
difference  cannot  be  accounted  for  by  any  reasonable  adjustment  for 
poor  mixing.  It  is  possible  that  the  faster  rate  observed  is  the 

2 

result  of  the  quenching  of  some  precursor  which  is  the  dominant  N(  D) 
source  in  their  higher  pressure  system  (~  30  Torr) .  The  non-linearity 
observed  in  some  of  their  plots  would  tend  to  support  this 
hypothesis . 
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3 

The  method  of  generating  0(  P)  atoms  necessarily  means  that  NO 

molecules  will  be  present  in  the  flow  for  some  finite  time.  As 
kNC)turns  out  to  be  nearly  two  orders  of  magnitude  larger  than 

kQ,  this  will  always  be  significant  although,  by  carefully 


controlling  the  experimental  conditions,  the  NO  contribution  can 
be  kept  to  a  minimum.  It  is  essential  that  che  flow  conditions  be 
precisely  characterised  for  each  kinetic  run,  and  that  the 
analysis  carried  out  with  extreme  care,  if  meaningful  results  are 
to  be  recovered.  This  also  applies  to  the  NO  results.  In  this 
instance,  the  requirement  that  the  measurements  be  carried  out  at 
very  short  reaction  times  brings  us  to  the  limits  of  flow  tube 
technology.  The  mixing  time,  using  standard  injectors  and  flows, 
becomes  a  significant  fraction  of  the  total  flow  time  and  must  be 

carefully  considered  to  avoid  reporting  artificially  low  values. 
The  accuracy  of  the  reported  value  of  kNQ  is  primarily  dictated 

by  the  confidence  limits  associated  with  k.,.  Should  a  more 

precise  value  of  k^  become  available,  the  confidence  limits  on 

kNQ  may  be  reduced  accordingly. 

Experimental  measurements  of  kQ  encounter  a  number  of 


difficulties,  the  most  immediate  of  which  are,  the  generation  of 

3 

a  pure  source  of  0(  ?)  atoms,  their  transportation  to  the 
reaction  zone  and  the  subsequent  efficient  mixing  of  the 
reactants.  The  method  of  titrating  NO  directly  into  the  flow  tube 

ingeniously  solves  the  transportation  and  mixing  difficulties  but 

,  3 

introduces  an  additional  complication  as  the  0(~P)  precursor 


itself  is  a  highly  reactive  impurity.  With  precise 
characterisation  of  flow  conditions  and  careful  analysis  of 
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results,  an  accurate  value  for  can  be  recovered. 

Recent  models  for  atmospheric  NO  concentrations  have  tended  to 
favor  larger  values  of  k^,  and  it  would  appear  that  experiment  and 

theory  have  finally  reached  a  consensus  regarding  the  preferred  value 

of  this  rate  coefficient.  Bates34  has  shown  that  the  fac> :  rate 

coefficient  reported  by  Jusinski  et  al.11  is  incompatible  with 

observed  data  on  01  red  line  emission  at  630  nm.  Also,  the  most 

35 

recent  modeling  calculation  by  Fesen  successfully  accounts  for 

-12  3  -1 

observed  NO  distributions  using  a  value  of  kQ  =  1x10  cm  s  ,  which 
coincides  with  the  most  recent  experimental  determinations. 
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Table  I.  Kinetic  Data  for  Determination  of  N(2D)  +  NO  -»N2+  0  Rate 
Coefficient 


t,  ms 

lower  ] 

bound 

knof  cm3  molecule-1  sec-1 
mid-range 

x  10-11 

upper 

bound 

1.82 

6.11 

± 

0.18 

7.72 

± 

0.22 

8.79 

+ 

0.26 

6.65 

± 

0.28 

8.52 

± 

0.36 

9.73 

+ 

0.41 

2.625 

7.38 

± 

0.31 

8.77 

± 

0.37 

10.30 

+ 

0.43 

7.53 

+ 

0.30 

8.88 

± 

0.35 

10.60 

± 

0.42 

3.175 

5.81 

± 

0.22 

7.13 

± 

0.27 

8.89 

+ 

0.34 

6.64 

± 

0.40 

8.24 

+ 

0.50 

10.40 

± 

0.63 

3.385 

6.56 

0.30 

7.84 

± 

0.36 

9.47 

+ 

0.43 

6.56 

+ 

0.26 

7.76 

+ 

0.31 

9.28 

+ 

0.37 

3.865 

4.39 

+ 

0.18 

5.13 

± 

0.22 

6.06 

+ 

0.25 

4.61 

+ 

0.24 

5.38 

± 

0.28 

6.38 

+ 

0.23 

3.895 

4.61 

± 

0.13 

5.67 

± 

0.16 

7.15 

+ 

0.20 

Weighted 

5.00 

5.49 

± 

± 

0.18 

0.07 

6.13 

6.70 

± 

± 

0.22 

0.08 

7.63 

8.06 

+ 

± 

0.27 

0.10 
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Table  II.  Kinetic  Data  for  Determination  of  N(2D)  +  0(3P)  Rate 
Coefficient 


t ,  ms 

lower 

ko ,  cm3 
bound* 

molecule-1  sec-1  x  10-13 
mid-rangeb 

upper 

boundc 

23 

4.74 

+ 

0.37 

7.60 

i 

0.37 

9.17 

+ 

0.37 

3.37 

+ 

0.50 

6.38 

± 

0.49 

8.03 

+ 

0.49 

31.7 

5.33 

+ 

0.23 

6.81 

± 

0.23 

7.57 

± 

0.23 

5.08 

+ 

0.31 

6.78 

i 

0.31 

7.64 

• 

+ 

0.31 

36.1 

7.12 

+ 

0.26 

8.51 

± 

0.26 

9.27 

+ 

0.26 

6.01 

± 

0.24 

7.44 

i 

0.24 

8.27 

+ 

0.24 

41.4 

6.62 

± 

0.21 

7.58 

± 

0.21 

8.14 

± 

0.21 

6.01 

± 

0.26 

6.99 

± 

0.26 

7.53 

+ 

0.26 

CO 

5.46 

± 

0.28 

6.21 

± 

0.27 

6.63 

+ 

0.27 

5.32 

+ 

0.28 

6.09 

+ 

0.28 

6.52 

+ 

0.28 

55 

6.19 

+ 

0.17 

6.65 

± 

0.17 

6.90 

+ 

0.17 

5.54 

+ 

0.23 

6.05 

± 

0.23 

6.33 

+ 

0.23 

Weighted 

mean 

5.85 

± 

0.08 

6.93 

± 

0.08 

7.64 

± 

0.08 

aUsing  kjj 
bUsing  kn 
cUsing  k^ 


2.5xl0-11 
3 . 4xl0-11 
4 . 3xl0-11 


cm3sec-1, 

cm3sec-1 

cm3sec-1 


kjvjo  =  5.49xl0-11  cm3sec_1 
kNo  =  6 . 70x10"^  cm3sec_1 
kfjo  =  8.06xl0-11  cm3sec-l 
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Table  III.  Experimental  and  a  Representative  Sample  of  Modeling 
Determinations  of  k0 


drafts."1 

Experimental  (1.8±0.6)xl0  -12 

<1.8  x  1CT12 
(2 .1±0 .8)  xlO-11 
(1.06±0.26)xl0“12 
(6. 9+1.0)  xlO"13 
Modeling  l.xlO-12 

<10'14 

4.10- 13 
l.xlO'12 
4.10"13 

6 . 10- 13 
~lxl0-12 


Source 

Reference 

Davenport  et  al. 
(1976) 

33 

Ianuzzi  & 

Kaufman  (1980) 

9 

Jusinski  et  al . 
(1988) 

11 

Piper  (1988) 

32 

this  work 

Strobel  gt  ,al. 
(1970) 

2 

Oran  gJL_al- 
(1975) 

12 

Frederick  & 

Rusch  (1977) 

13 

Cravens  et  al. 
(1979) 

14 

Rusch  &  Gerard 
(1980) 

15 

Richards  gt— al- 
(1981) 

36 

Fesen  gt_al- 

35 

(1989) 


Reaction  Time/ms 


Figure  2.  Variation  in  the  effective  rate  coefficient, 

(5))  with  reaction  time.  The  reaction  time  was  varied  by 
adjustments  to  the  injector  position  and  the  flow  speed. 

4 

Consequently,  the  flow  conditions  (i.e.  total  pressure,  [N(  S) ] 
etc.)  varied  with  the  reaction  time. 


Figure  3.  Variation  in  the  fractional  amount  of  initially  added 

that  remained  unreacted  at  the  MPI-detection  point.  As  with  the 
,  measurements,  flow  conditions  varied  with  reaction  time. 


F 1  a wt i me /ms 


2 

Figure  4.  Variation  in  the  decay  coefficient  (  -dln[N(  D) ] 

2 

/d[02])  with  reaction  time  for  N{  D)  quenching  by  €>2  •  Flowtube 

pressure  was  0.46  Torr  at  a  linear  flow  velocity  of  6700  cm  s-1-  ' 
error  bars  represent  the  statistical  limits  to  the  decay 
coefficient  obtained  from  a  weighted  linear  regression  fit  to  the 
individual  data  sets. 


[NO]dt,  which  was  calculated  by  numerical  integration  of  Eq 
.t  =  3.90ms  (corrected  for  0.48ms  mixing  time),  kM  =  3.4x10 
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Figure  7.  Decay  of  the  N(2D)  MPI-signal  as  a  function  of 
J0  !0]dt.  The  bottom  trace  was  calculated  assuming  instantaneous 
conversion  of  [NO]  to  [0]  U.e.,  £  £0J dt- |N0)ot)  and  therefore  no 

correction  for  depletion  of  N(2D)  by  NO  is  necessary.  A  linear 
regression  fit  to  these  data  gives  a  value  of  kn  =  2.03xl0-*2 

cm  s  ,  The  upper  trace  was  calculated  using  Eq.  (15)  with  kN  * 
3.4x10  11  cm3s  \  kNQ  =  6.7xlO"U  cm3s~1,  [N(4S)J  =  3.28mTorr  and  t 
23.0ms.  This  gives  a  value  of  kQ  =  7.60X10~13  cmV1.  The  difference 

in  the  two  plots  shows  the  relative  magnitude  of  the  NO  contribution 
to  the  N(2D)  signal  decay. 
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3.  Reactions  of  N(2P) 

3.1.  N(2P)  Kinetics 

2 

N(  P)  kinetics  may  be  investigated  using  an  identical  experimental  set-up  to  that  used  for 

2  2  1 
N(  D)  studies.  MPI  detection  of  N(  P)  utilizes  the  two  photon  transition 

N[3p  2P°j]  <-  N[2p3  2P°j]  (3.1) 

2 

at  290  nm  with  a  third  photon  then  ionizing  the  atom.  Measurements  of  rates  for  the  reaction  of  N(  P) 

3 

with  C>2  and  NO,  and  for  quenching  by  0(  P)  were  taken.  The  experimental  procedure  and  kinetic 

2 

analysis  are  identical  to  those  given  in  the  previous  section  for  the  N(  D)  measurements. 

3.2.  Results 

3.2.1.  N(2P)  +  02  Experiments 

2 

A  previous  measurement  of  the  rate  coefficient  for  the  reaction 

N(2P)  +  02 ->N0  +  0  (3.2) 

-12  3  -1 

made  using  the  FACELIF  apparatus,  reported  a  value  of  1.8x10  cm  s  .  This  was  somewhat 
lower  than  some  previously  reported  values'3’  .  We  therefore  undertook  a  remeasurement  of  this  rate 
coefficient  in  which  the  effects  of  the  finite  mixing  time  were  explicitly  accounted  for.  This  was  done, 
as  reported  in  the  previous  section,  by  performing  rate  measurments  over  a  range  of  reaction  lengths, 
under  uniform  flow  conditions.  These  measurements  were  carried  out  under  fast  flow  conditions  and 


the  quenching  rate  for  each  data  set  plotted  as  a  function  of  the  flow  time.  The  results  are  plotted  in 

12  3  1 

Fig.  3.1.  The  slope  of  this  graph  yields  a  value  for  the  rate  coefficient  of  3.4x10  cm  s  .  The  x- 
axis  intercept  gives  a  value  for  the  effective  mixing  time  for  the  particular  flow  conditions  which  in 


this  instance  is  approximately  1.2  ms.  Neglect  of  this  effect  leads  to  a  considerable  underestimation  of 

flip  fptp  POPfflPIPnt  ' I  Vl A  rPTV„*tP/1  xioImo  1C  in  (rn/vl  nnroftmont  iiii.lx  nrainfinii  *-n n o f. . . »H n .1  * (1 3 ?4 
v»»v  *uw  vwuiivtvint  Aits,  tuiuv  xo  ah  gwu  u^ivuiiivni  mui  pivriuuo  uiuuouivuiuud  ♦ 
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3.2,2.  N(2P)  +  NO/O 

As  with  N(2D),  the  variation  in  the  effective  rate  coefficient  (the  rate  coefficient  calculated 

3 

assuming  that  all  the  initially  added  NO  is  converted  instantaneously  to  0(  P) )  with  reaction  time  was 

measured  over  a  wide  range  of  flow  times.  The  results  are  shown  in  Fig.  3.2.  The  error  bars  on  each 

data  point  represent  the  statistical  errors  obtained  from  a  linear  regression  fit  to  the  individual  data 

2 

sets.  The  temporal  behavior  is  identical  to  that  observed  for  N(  D)  (see  Fig.  2  in  Section  2).  The 
reaction 

N(2P)  +  NO  ->  N2  +  0  (3.3) 

5  2 

is  known  to  be  rapid  .  Therefore,  as  for  N(  D),  the  short  reaction  time  measurements  may  be 
considered  to  give  a  lower  limit  on  the  rate  coefficient  for  reaction  (3.3).  It  might  also  be  expected  that 
the  longer  reaction  time  measurements  will  give  an  upper  limit  to  the  N(2P)+0  quenching  rate 
coefficient. 

2  3 

In  order  to  derive  rate  coefficients  for  N(  P)  reaction  with  NO  and  quenching  by  0(  P),  it  is 
necessary  to  know  the  N(4S)  atom  concentration.  For  the  NO  measurements,  [N(^S)]  was 
determined  using  a  variation  of  the  technique  described  in  the  previous  section.  The  equation 
describing  the  NO  decay  along  the  flowtube  is 

[NO]  =  [NO]0exp(-kN[N]t)  (3.4) 

provided  that  [NO]  « [N].  In  this  instance,  the  NO  7(0,0)  MPI  signal  was  monitored  as  a  function 

of  the  injector  position.  This  is  equivalent  to  varying  t  and  therefore  a  plot  of  In  [MPI-signal]  vs.  t  will 
have  a  slope  of  kj^[N].  k^  is  known  from  the  literature^  and  thus  [N]  is  derived  for  all  points  along 

the  flowtube.  The  results  are  plotted  in  Fig.  3.3. 

The  results  do  not  display  the  linear  behavior  that  might  be  expected  over  the  whole  range. 
Three  distinct  regions  are  apparent.  The  two  shorter  reaction  time  regions  (t  <  3.2ms)  indicate  that  [N] 
is  different  for  these  two  ranges  of  injector  position.  [N]  measurements  made  using  the  technique 
described  in  the  previous  section  and  direct  MPI  measurements  (next  section)  confirm  that  this  is  the 
case.  The  variation  in  [N]  would  appear  to  be  the  result  of  a  physical  effect  with  the  loop  injector 
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disturbing  the  spatial  distribution  of  the  afterglow  as  it  is  withdrawn  towards  the  point  at  which  the 

flowtube  diameter  increases  (see  Fig.  1  in  Section  2).  While  the  effect  appears  to  have  a  sharp  onset  it 

is  probable  that  there  is  some  variation  in  [N]  over  the  whole  range.  The  main  source  of  error, 
however,  is  expected  to  be  from  the  error  bars  associated  with  the  literature  value  of  kj^.  The  values 

obtained  are,  for  t  <  2.7  ms,  [N]  =  (0.21 1  +  0.055)  mTorr  and,  for  3.2  <  t  <  2.7  ms,  [N]  =  (0.41 1 
&  0.107)  mTorr. 

As  the  injector  is  drawn  further  back  towards  the  expansion  point,  the  flow  is  disturbed  such 
that  the  main  and  quencher  flows  do  not  become  fully  mixed  for  some  distance  beyond  the  addition 
point.  The  net  result  is  that  the  effective  contact  time  between  the  reactants  remains  approximately 
constant  over  the  whole  range.  Care  must  be  taken  to  avoid  this  region  when  taking  kinetic 
measurements. 

The  results  of  three  kinetic  runs  were  computed  using  the  analysis  given  in  the  previous 

2 

section.  The  effective  rate  coefficient  measurements  (Fig.  3.2)  suggest  that,  as  for  N(  D),  the  O-atom 

contribution  to  the  decay  may  be  neglected  for  fast  flows.  The  flow  conditions  used  were  identical  to 

those  used  for  N(^D)  +  NO  experiments  and  therefore  the  same  effective  mixing  time  of  (0.48  ± 

0.19)  ms  was  used.  The  data  are  presented  in  Figs.  3.4-3.6.  The  rate  coefficients  obtained  using  this 

procedure  are  (4.64  ±  0.12)xl0_1 1,  (4.52  ±  0.26)xl0_1 1,  and  (4.32  ±  0.12)xl0_1 1  cm3s'1  for  t  = 

1.78,  1.98,  and  2.16  ms,  respectively. 

2 

This  value  for  the  N(  P)  +  NO  reaction  rate  coefficient  was  then  used  to  calculate  a  value  for 

the  O  atom  quenching  rate  coefficient.  As  for  N(^D),  these  measurements  are  best  taken  at  long 

flowtimes  with  high  [N].  The  results  of  four  measurements  are  shown  in  Figs.  3.7-3.10.  These  yield 

rate  coefficient  values  of  (2.45  ±  0.07)xl0‘12,  (3.05  ±  0.13)xl0'12,  (2,83  ±  0.04)xl0'12,and  (4.29 
-12  3  -1 

±  0.14)xl0  cm  s  ,  for  t  =  38.1, 31.5, 27.7,  and  22.1  ms,  respectively.  A  number  of  other 
measurements  were  also  taken.  These  all  exhibited  a  trend  towards  lower  rate  coefficient  values  as  the 


reaction  time  increased. 


42 


3.3  Discussion 

It  has  recently  been  reported^  that  microwave  discharges  do  not  provide  a  clean  source  of 
N(  P).  The  results  reported  by  Piper  exhibit  a  rapid  initial  decay  followed  by  a  considerably  slower 
decay  for  larger  quencher  additions.  This  was  interpreted  as  the  initial  rapid  quenching  of  some 

unidentified  precursor  to  N(  P),  followed  by  quenching  of  the  remaining  microwave  generated 

2  2 
N(  P).  It  was  suggested  that  this  second  region  represented  the  true  N(  P)  rates. 

The  longer  reaction  time  measurements  showed  indications  of  non-exponential  behavior 
in  some  instances  (see  Figs.  3.7  -  3.10).  Also,  the  recovered  rate  coefficients  decreased  with 
increased  reaction  time  even  after  correction  for  the  NO  contribution.  The  range  of  values  indicated  by 
these  results  is  considerably  lower  than  previously  reported  O-atom  quenching  rate  coefficients, 
which  range  from  1  x  10'11  to  7  x  10'11  cm3  s'1 7>8.  These  results  clearly  indicate  that,  for  these  time 
scales,  we  are  net  observing  a  pure  exponential  decay  of  discharge  generated  N(2p). 

3.4.  Conclusions 

2 

It  is  not  possible  to  infer  anything  regarding  the  N(  P)  +  O  quenching  rate  coefficient  from 

these  results.  The  variation  in  the  effective  rate  coefficient  (Fig.  3.2)  with  time  is  clearly  misleading. 

2  6 
The  N(ZP)  +  NO  results  were  therefore  recalculated  using  the  most  recently  reported  value  for  the 

2  -113-1 

N(  P) +  0  quenching  rate  coefficient  of  1.7x10  cm  s  .  The  data  appear  to  give  reasonable 

exponential  decays  and  consistent  rate  coefficient  values.  Given  the  short  reaction  time  and  low 

[N(^S)],  the  O-atom  correction  is  still  comparatively  minor.  This  is  illustrated  graphically  in  Fig.  3.1 1 

where  the  magnitude  of  the  correction  is  on  the  order  of  5%.  The  recalculated  values  are  (4.38  ± 

0.12)xl0"11,  (4.24±0.25)xl0_11  and  (4.00  ±0.12)xl0‘U  cmV1  for  t  =  1.78,  1.98,  and  2.16 

ms,  respectively,  giving  a  mean  value  of  (4.21  ±  0.17)xl0  cm  s  .The  quoted  limits  here  only 

include  errors  from  the  fitted  slopes  and  take  no  account  of  uncertainties  in  [N]  or  the  mixng  time. 

3  9 

This  value  is  approximately  25%  higher  than  previously  reported  values  ’  .  Given  the  limited 
dynamic  range  of  these  experiments,  it  is  not  possible  to  comment  on  the  suggestion^  that  this 
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2 

measurement  in  fact  reflects  the  reaction  rate  of  NO  or  O2  with  an  unidentified  N(  P)  precursor. 

Using  the  new  FACELIF  configuration  (see  Section  7),  incorporating  a  movable  MPI  detector,  it 

should  be  possible  to  observe  the  N(  P)  decay  along  the  flowtube.  Any  deviation  from  the  predicted 

2  6 

exponential  time  decay  (governed  by  wall  diffusion)  will  indicate  the  presence  of  N(ZP)  precursors0 

2 

or  a  possible  N  atom  recombination  source  . 
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figure  3.2.  Variation  in  the  quenching  coefficient  (-d  ln[N(2P)]/d[C>2])  with  flowrime.  Flowtube 
pressure  =  0.39  Torr,  flow  velocity  =  7400  cm  s'1. 
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Figure  3.3.  Variation  in  [NO]  with  injector  position  (=  flowtime).  Flowtube  pressure  =  0.38Torr, 
flow  velocity  =  7400  cm  s'*  and  [NO]n  =  0.025  mT. 


48 


(B[(d_)N]/C(d, )N3  )ul 
z  z 


Figure  3.5.  As  Fig.  3.4  with  t  =  1 .98ms,  [N]  =  0.21 1  mT  and  kNO  =  (4.52  ±  0.26)xl0' 
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(B[(d2)N]/C(d2)N]  )ul 


Figure  3.6.  As  Fig.  3.4  with  t  =  2.16  ms,  [N]  =  0.21 1  mT  and  kNn  =  (4.32  ±  0.16)xl0' 


so 


-signal  as  a  function  of  /[Ojdt.  Flowtube  pressure  =  4.1  Torr, 
=  (4.29  ±  0.14)xl0"12  cm3  s'1. 


pressure  =  3.0  Torr,  [N]  =  3.06  mT,  t  =  27.7  ms  and  k( 
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Figure  3.9.  As  Fig.  3.7  with  [NT  =  2.87  mT,  t  =  31.5  ms  and  kQ  =  (3.05  ±  0.13)xl0 
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(0C(d2)N3/C(d2)N3  >ui 


Figure  3.10.  As  Fig.  3.7  with  [N]  =  3.54  mT,  t  =  38.1  ms  and  kQ  =  (2.45  ±  0.07)xl0 
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4.  RESEARCH  NOTE  [Spectrochiniica  Acta,  v.  46A,  1990] 


DETECTION  OF  N(4S)  BY  RESONANTLY  ENHANCED  MULTIPHOTON  IONIZATION  SPECTROSCOPY 


(Received  1  December  1989) 


Abstract  -  We  report  the  First  instance  of  detection  of  N(4S)  atoms  by  Resonantly  Enhanced 
Multiphoton  Ionization  (REMPI).  The  magnitude  of  the  ion  signal,  after  correction  for  laser  power 
dependence,  is  directly  proportional  to  the  N(4S)  concentration  and,  once  calibrated,  is  a  direct  and  highly 
sensitive  method  for  N(4S)  concentration  measurements  in  fast-flow  kinetics  experiments. 


INTRODUCTION 

In  experiments  utilizing  active  nitrogen 
flows,  it  is  generally  necessary  to  know  the 
ground  state  nitrogen  atom  concentration.  This 


value  is  usually  determined  using  a  standard  NO 
titration  technique  (1).  This  method  makes  use  of 

the  reactions 

N(4S)  +  N0~»N2  +  0(3P) 

(1) 

and 

NO  +  O  +  M  ->  N02*  +  M 

(2) 

Reaction  (1)  is  rapid  (2)  and  serves  to  deplete 
(N(4S)]  stoichiomctrically.  Monitoring  the  N2 
first  positive  emission  around  580nm  gives  a 
fluorescence  signal  proportional  to  (N(4S)]2  for 
[NO]<[N(4S)].  Beyond  the  end  point  N02* 
emission  following  reaction  (2)  gives  a  signal 
proportional  to  NO  in  the  same  spectral  range. 
This  technique  works  well  provided  that  reaction 
(1)  has  sufficient  time  to  go  to  completion  and 
N(4S)  is  present  in  sufficient  quantity.  If  this  is 
not  the  case  then  both  the  N2  first  positive  and 
N02*  emission  will  contribute  to  the  signal 
level. 

Many  reactions  involving  mctastable 
atoms  and  molecules  require  fast  flow  conditions 
due  to  the  rapid  decay  of  the  metastable 
concentration.  The  reactions 


N(2D)+N0-»N2  +  0 

(3) 

N(2P)+02— >NO  +  0 

(4) 

are  examples  of  this  [3-5].  In  order  to  calculate 
the  rate  coefficient  of  reaction  (3),  (N(4S)]  must 
be  known  in  order  to  determine  the  rate  of 
depletion  of  NO  via  reaction  (1).  Under  certain 
conditions,  the  flow  time  and  N(4S)  concentration 
may  be  such  that  reaction  (1)  docs  not  go  to 
completion  and  the  standard  NO  titration 
technique  will  no  longer  give  accurate  values  for 
[N(4S)j,  The  use  of  Resonantly  Enhanced 
Multiphoton  Ionization  (REMPI)  to  measure 
relative  N(^D)  and  N(^P)  concentrations  is  well 
established  (5,6).  Here  we  report,  to  our 
knowledge,  the  first  instance  of  REMPI  detection 


of  N(^S)  atoms.  The  magnitude  of  lire  ion  signal 
is  directly  proportional  to  the  N(**S) 
concentration  and,  once  calibrated,  may  be  used 
to  measure  N(^S)  absolutely. 

EXPERIMENTAL 


These  investigations  were  earned  out 
using  the  FACELIF  afterglow  reactor  which  has 
been  described  in  detail  elsewhere  [4,7].  A 
schematic  of  the  experimental  apparatus  is  shown 
in  Fig.  1.  The  N(4S)  atoms  are  generated  in  a  flow 
of  1-2%  N2  in  He,  using  a  microwave  discharge. 
The  current  created  by  the  photoions,  across  a 
pair  of  Ni  wire  grids,  biased  at  90V,  is  amplified 
using  an  Ithaco  model  121 1  cuncm  amplifier  and 
detected  on  an  EG  &  G  model  162  Boxcar 
Averager.  The  laser  power  is  also  monitored 
using  a  Molcctron  Pyroelectric  joulcmctcr  (J3- 
09).  The  fluorescence  measurements  arc  made 
using  a  photomultiplier  tube  through  a  580  ±  5 
nm  filter.  The  NO  used  for  these  titrations  is  a  5% 
mixture  of  NO  in  He  and  is  added  through  a  loop 
injector  upstream  of  the  detection  region. 
Typical  flow  times  were  25-30ms  with  total  lube 
pressures  of  2.5-3.5  Torr.  These  slow  flow/high 
pressure  conditions  arc  best  suited  to  making  NO 
titration  measurements  of  [N(4S)]  giving 
comparatively  intense  signals  and  allowing 
reaction  (1)  ample  lime  to  go  to  completion.  The 
N(4S)  concentration  could  be  varied  cither  by 
adjusting  the  N2  to  He  ratio  (this  was  typically  1'- 
3%  N2  "in  He)  or  more  conveniently  by  altering 
the  microwave  power. 

The  transition  N(2p23p  4Dj°)  <—  N(2p^ 
4S3/2°)  may  be  resonantly  excited  by  two  21 1  nm 


photons  (8),  and  the  atom  then  ionized  by  a  third 
photon  of  the  same  energy.  Laser  radiation  at 
211  nm  may  be  generated  by  stimulated  Raman 
frequency  shifting  of  the  doubled  output  of  a 


gives  pulse  energies  of  less  than  100  mJ. 
Significantly  higher  pulse  energies  arc  desirable 
for  REMPI  detection.  Here  we  have  used  a  BBO 
crystal  to  double  the  output  from  a  Nd:YAG 
pumped  dye  laser  (Quanta-Ray  DGR2A,  PDL-1, 
INRAD  Auioiracker  II)  directly.  Stilbcnc  420 
laser  dye  gave  pulse  energies  of  approximately 
0.5  mJ,  measured  after  exiting  the  flow  tube. 
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This  power  level  decreased  rapidly  in  time  with 
the  dye  having  a  useful  lifetime  of  about  20 
minutes.  The  decrease  proved  advantageous  when 
taking  measurements  of  the  power  dependence  of 
the  ion  signals,  enabling  a  large  number  of  data 
points  to  be  taken  without  the  extensive  use  of 
neutral  density  niters  to  attenuate  the  beam. 
When  making  quantitative  measurements  of 
[N(4S)j ,  however,  it  is  desirable  to  keep  the  laser 
power  as  constant  as  possible  to  minimize 
corrections  necessary  for  power  normalization  of 
ion  signals.  We  found  that  the  use  of  ~  four  times 
larger  dye  reservoirs  greatly  enhanced  the  useful 
life  of  the  dye  and  also  minimized  power 
variations  during  the  course  of  data  acquisition. 

ABSOLUTE  [N(4S)]  MEASUREMENTS 

The  N(4Dj°)  state  is  split  into  four 
levels,  J=  1/2,  3/2,  5/2,  and  7/2,  at  94770.85, 
94793.46.  94830.86,  and  94881.79  cm'1, 
respectively.  The  relative  intensities  are  simply 
given  by  the  upper  slate  degeneracy  (2J+1)  and 
the  observed  spectrum  will  be  a  quartet  with 
intensity  ratios  of  1:2:3:4  (J=l/2:3/2:5/2:7/2). 
The  observed  spectrum  is  shown  in  Fig.  2.  The 
most  intense  line  was  used  throughout  for  all 
quantitative  measurements  of  [N(4S)1.  To 
normalize  the  signal  levels  accurately  for 
variations  in  laser  power,  it  was  necessary  to 
characterize  the  power  dependence  of  the  ion 
signal  precisely  prior  to  any  data  acquisition. 
The  ion  signal  is  proportional  to  some  power  N 
of  the  laser  intensity  /,  and  therefore  a  plot  of  In 
(ion  signal)  vs.  In  /  will  have  a  slope  equal  to  N. 
The  results  of  doing  this  are  shown  in  Fig.  3.  The 
slope  varies  between  2.8  and  1.5  going  from  low 
to  higher  power.  This  indicates  that  the  two 
photon  transition  is  becoming  saturated  at  the 
higher  powers.  A  similar  power  dependence  was 
observed  for  (2+1)  photon  ionization  of  N(2D)  at 
269nm  [10],  although  somewhat  different 
conclusions  were  drawn  in  that  work  regarding 
the  relative  saturation  of  the  excitation  and 
ionization  steps.  The  two  photon  excitation 
cross-section  for  this  transition  has  been 
calculated  (9]  and  therefore  the  transition  rate 
may  be  calculated.  The  fit  of  the  data  in  Fig.  3  is 
obtained  using  a  simple  rate  equation  model  for  a 
three  level  system  (11).  Strictly  speaking, 
correct  normalization  of  signals  for  power 
variations  requires  that  signal  and  laser  power  be 
recorded  for  each  pulse  (12).  Here  we  have 
recorded  the  average  ion  signal  and  laser  power 
over  a  large  number  of  pulses  with  the  beam 
tightly  focused  between  the  detector  grids.  The 
resulting  ion  signal  therefore  consists  of 
contributions  over  a  range  of  saturation 
conditions.  It  is  thus  inappropriate  to  comment 
on  the  magnitudes  of  the  experimentally 
observed  rates  for  either  the  ionization  or  the  two 
photon  excitation.  For  the  purposes  of  this  study 
it  is  only  necessary  to  be  able  to  understand  the 
bulk  behavior  of  the  system.  Fig.  3  indicates 
that  for  small  variations  in  average  laser  power 
(±5%)  it  is  a  reasonable  approximation  to  use  the 


same  laser  power  dependence  factor  when 
normalizing  the  ion  signals. 

For  conditions  where  reaction  (1)  has 
time  to  go  to  completion,  (N(4S))  may  be 
determined  by  observing  the  reduction  in  the  ion 
signal  as  NO  is  added.  The  x-axis  intercept  on  a 
plot  of  the  ion  signal  vs.  [NO)  gives  the  value  of 
[N(4S)J.  The  slope  of  this  line  serves  to  calibrate 
the  ion  signal  level  in  terms  of  absolute  [N(4S)| 
concentrations.  In  order  to  verify  the  accuracy  of 
this  technique  for  [N(4 S ) ]  measurement, 
comparative  measurements  were  taken  using 
REMPI  and  the  standard  fluorescence  monitoring 
titration  technique  [1].  The  results  are  listed  in 
Table  1.  The  REMPI  results  are  the  x-axis 
intercept  of  a  weighted  linear  regression  fit  of  NO 
titration  data  sets  with  ±2o  error  limits.  The 
standard  titration  measurements  may  be 
considered  accurate  to  ±5%.  The  results  arc  in 
good  agreement  indicating  that  the  REMPI 
technique  yields  accurate  values  for  [N(4S)). 

DISCUSSION 

REMPI  provides  a  sensitive  method  for 
quantitatively  detecting  ground  state  nitrogen 
atoms.  Number  densities  on  the  order  of  10'° 
atoms  cm'3  were  easily  detected  compared  to  a 
limit,  for  our  flow  system,  of  around  2xl013 
atoms  cm'3  using  the  more  usual  titration 
technique.  Below  this  level  the  fluorescence 
measurements  encountered  problems  in  the  form 
of  stray  light  from  the  discharge  reaching  the 
detector.  A  more  serious  problem  can  arise  from 
fluorescence  due  to  sources  other  than  N(4S)  atom 
recombination,  which  contributes  significantly 
to  the  signal  intensity  [4,13].  Unlike  the 
standard  technique,  REMPI  offers  a  state  specific 
method  of  monitoring  [N(4S)]  directly  and  is 
therefore  less  susceptible  to  extraneous  signals. 
The  signal  level  may  be  calibrated  under  slow 
flow/high  [N(4S)J  conditions.  The  magnitude  of 
the  ion  signal  may  then  be  used  to  give  a  direct 
measurement  of  [N(4S)1  for  experiments  requiring 
fast  flows  where  [N(4S)|  may  be  low;  these 
condition  are  inappropriate  for  standard  titration 
measurements. 
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Table  1.  Comparison  of  (N(4S)]  titrations  by  NO  using  standard  chemiluminescence  detection  and  direct 
MPI  detection  under  flow  conditions.  (N(4S)]  given  in  mTorr. 


Standard  NO  Titration 


MPl-TItratlon 


4.38 

3.28 

2.65 

2.02 

1.97 

1.42 


4.48  ±0.18 
3.40  +  0.10 

2.50  ±0.14 

1.96  ±0.12 

1.97  ±0.12 

1.51  ±0.07 


FIG.  1.  Schematic  of  FACELIF  experimental  apparatus.  The  solid  circle  represents  a  right  angle 
connection  to  tire  pumping  system. 
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5  «A  Kinetic  Model  for  Resonantly  Enhanced  Multiphoton  Ionization  of  NO 

I.  INTRODUCTION 

Multiphoton  ionization  (MPI)  is  one  of  the  most  sensitive  techniques  for  state  specific  detection 
of  gas  phase  molecules.  For  accurate  treatment  of  the  population  distributions  measured  by  MPI.  the 
dependence  of  the  ion  signal  on  laser  intensity  and  the  state  specific  transition  must  be  known.  The 
dynamics  of  MPI  are  complicated.  Formal  density  matrix  treatment  and  inclusion  of  all  the  nonlinear 
processes  would  be  cumbersome,  to  say  the  least.  Fortunately,  for  most  experimental  conditions,  a 
more  intuitive  rate  equation  approach  may  be  taken.’*7  Such  is  the  case  for  (1  +  1)  MPI  of  NO. 

The  availability  of  tunable  u.v.  lasers  at  226  nm  and  the  low  ionization  potential  of  NO  make 

the  NO  y-band  2 E + <-2II  an  ideal  system  for  (1  +  1)  MPI  spectroscopy.6*12,14  From  a  rate  equation 

model  of  MPI  .  a  two  photon  process  would  depend  quadratically  on  the  intensity  of  radiation.  But. 
for  a  sufficiently  intense  source,  the  ground  to  intermediate  state  resonance  will  saturate,  ionization 
will  be  the  rate  limiting  step,  and  the  dependence  will  be  linear.11  Because  of  this  linearity,  many  in¬ 
vestigators  try  to  work  in  this  regime.  Shot  to  shot  variations  in  laser  intensity  may  be  corrected  by 
averaging  the  laser  power  via  boxcar.  But  this  linear  region  is  really  a  transition  between  quadratic 
and  zeroth  order  behavior  (see  Section  III  and  ref.  2).  The  laser  power  necessary  for  first  order 
dependence  will  vary  for  different  vibrational  subbands  and  will  even  vary  among  rotational 
transitions.14 

Zare  and  Jacobs'^'14  have  proposed  an  excellent  method  for  empirically  fitting  the  ion  signal  to 
integrated  laser  intensity.  Their  "isopower"  spectra  are  correctly  adjusted  for  variations  in  laser  inten¬ 
sity.  But  the  relationship  between  populations  of  different  vibrational  bands  is  still  ambiguous.  In 

this  paper,  we  will  show  a  simple  kinetic  model  based  on  the  spectroscopic  parameters  for  NO  that 

predicts  the  form  of  isopower  spectra. 

The  model  is  used  to  relate  ion  signals  from  different  vibrational  subbands  to  a  common  scale. 
We  will  use  these  results  to  discuss  the  correction  of  the  state  specific  kinetics  for  the 
N(4S)+O2+NO(v=0-7)  reaction  measured  via  MPI.12 

2.  METHODOLOGY 
2-1-THE  KINETIC  MODEL 

To  simulate  (1  +  1)  MPI  spectra,  it  is  necessary  to  define  the  frequency  at  which  transitions  oc¬ 
cur  and  the  magnitude  of  the  ion  signal.  The  resonant  ground  to  intermediate  state  energy  separation 
is  the  frequency  of  the  transition.  In  the  case  of  NO  A2E++X2II  these  frequencies  may  be  deter¬ 
mined  using  the  equations  given  by  Herzberg20  and  the  molecular  constants  for  NO.18  Of  more  inter¬ 
est  here  is  the  magnitude  of  the  ion  signal  and  its  dependence  on  laser  power.  The  ion  signal  is 
proportional  to  the  number  of  ions  produced  during  a  laser  pulse.  The  kinetic  model  assumes  a  tun¬ 
able.  pulsed  laser  produces  radiation  of  intensity  I  (mW/cm2  s*’)  and  pulse  duration  At  having  a  rec¬ 
tangular  temporal  width.  When  the  laser  is  tuned  to  a  specific  ro-vibrational  transition,  the  processes 
of  induced  absorption,  stimulated  emission,  and  ionization  will  proceed  for  the  duration  of  the  pulse. 
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By  solving  the  kinetic  equations  appropriate  to  the  transition  in  resonance,  the  fraction  of  the  NX' 
number  density  converted  to  ions  is  determined.  The  product  of  this  parameter  and  the  Boltzmann 
factor  for  the  transition  is  proportional  to  the  observed  ion  signal. 

It  is  instructive  to  consider  the  energy  level  diagram  for  the  NO  y-band  to  discern  the  nature 
of  the  transitions.  When  the  separation  between  2II3/2  and  2n,/2  is  large,  as  is  the  case  in  NO.  the 

multiplets  may  be  separated  into  two  subbands.  Selection  rules,  AJ  =  0,  ±1  and  allow  six 

transitions  in  each  subband,  twelve  total.  Transitions  which  originate  on  the  same  J"  but  terminate 
on  different  J’  within  the  same  K’  are  separated  only  by  the  spin  splitting  in  the  2E  +  .  This  separa¬ 
tion  is  very  small,  and  is  not  resolved  in  MPI.  For  the  kitu..v  model,  the  equations  for  the  four 
pair  of  coincident  transitions  must  be  solved  differently  than  for  the  noncoincident  transitions.  The 
equations  for  noncoincident  transitions  have  analytic  solutions  whereas  the  equations  for  coincident 
transitions  must  be  solved  numerically. 

For  the  case  of  noncoincident  transitions,  the  kinetics  are  described  by  a  three  level  system  (Fig 
5- la).  X  represents  the  population  in  a  J"  level  of  the  2II  state.  A  represents  the  population  in  a  J' 

level  of  the  2E+  state  corresponding  to  a  P,  Q.  or  R  transition.  C  represents  the  continuum  of  ion 

slates.  The  rate  constants  k)2  and  k21  are  for  induced  absorption  and  stimulated  emission.  k|0n  cor¬ 
responds  to  the  ionization  cross  section.  The  spontaneous  emission  rate.  I/T.  is  negligible  assuming 
a  5ns  laser  pulse,  but  is  included  for  completeness.  The  intermediate  state  lifetime  is  taken  as  216 
ns.21 

Equations  5-1  to  5-3  are  differential  equations  for  the  rate  of  change  of  population  in  X.  A. 
and  C. 


dX(t)/dt  =  -k,2X  +  k,,A  (5-1) 

dA(t)/dt  =  -(k|2  +  k.on  -  1/T)A  +  k,2X  (5-2) 

dC(t)/dt  =  ki0|,A  (5-3) 

Equations  5-1  and  5-2  are  simultaneously  solved  to  yield  expressions  for  X(t)  and  A(t). 

X(t)  =  (r2+k12/r,-r2)[er2'-ert']  (5-4) 

A(t)  =  {(r|  +  k12)(r2  +  k,2)/k21(r,-r2)}[er21-eril]  (5-5) 


r,  -  {-(m  +  k|2)  +  [(m  +  k|2)2-4(mk,2-k12k21)],/2}/2 
r2  =  {-(m  +  k12)-[(m  +  k12)M(mk12-k!2k2))]l/2}/2 


Equation  5-5  is  substituted  into  equation  5-3  for  a  solution  to  C(t) 

C(t)={k!ol)(r|  +  k.2)(r2+k12)/k2i(r1-r2)} 

{ (er2  - 1  )/r2)-(eri I )/r ,) }  (5-6) 

The  initial  conditions  set  A(0)  and  C(0)  equal  to  zero  and  X(0)  equal  to  one.  By  letting 
X(0)=l,  X(t),  A(t),  and  C(t)  at  all  later  times  are  the  fraction  of  the  NO  number  density  in  these 
levels. 
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For  the  case  of  coincident  transitions,  the  kinetic  system  includes  two  separate  intermediate 
levels  (Fig  5- lb).  The  differential  equations  for  this  four-level  system  are  Eqs.  5-7  to  5-10. 


dX(t)/dt= k2,  A+ k’2,  A’-(k,2+ k’  12)X  (5-7) 

dA(t)/dt = k,2X-(k2,  +  k|on  + 1  /  t)A  (5-8) 

dA'(t)/dt  =  k* ,  2X-(k’2,  +  k|0|(  + 1  /T)A’  (5-9) 

dC(.)/d.  =  kio|iA  +  k-.onA’  (5-10) 


An  analytic  solution  for  the  four-level  system  may  be  found  using  Laplace  transforms,  but  the 
equations  are  more  easily  integrated  numerically. 

To  complete  the  model  of  predicting  MPI  ion  signals  from  first  principles,  the  rate  coefficients, 
kp  and  k2),  must  be  determined  from  spectroscopic  parameters.  From  Einstein’s  theory  for  two-level 
systems  the  induced  absorption  and  stimulated  emission  rates  are  written: 

WI2  =  B12pwNj  (5-11) 

W2.  =  B2.P»N2  (5-12) 

W)2  and  are  the  rates  of  induced  absorption  and  spontaneous  emission.  B|2  and  B2|  are 
the  Einstein  B  coefficients  with  units  of  s'1  (Jem'"5 A'1)'1.  pw  is  the  energy  density  in  units  of 
Jcm'^A'1 .  The  rate  coefficients  k12  and  k2,  are  the  product  of  the  Einstein  B  coefficient  and  the 
energy  density. 

The  Einstein  B  coefficients  are  related  to  the  line  strength.  S.  of  a  ro-vibrational  transition  ac¬ 
cording  to  Eqs  5-13  to  5-14. 

B|2  =  (2n2X2/3h2e0cgj)S  (5-13) 

B2,  =  (2n2X2/3h2e0cg2)S  (5-14) 

g,  and  g2  are  the  electronic  statistical  weights  for  levels  1  and  2.  For  diatomics,  the  electronic 
statistical  weight  is  expressed  as 

(2-60^)(2S+l).  (5-15) 

The  Kronecker  delta,  &QA,  is  one  for  A  =  0.  that  is  Z  states,  and  zero  for  all  other  states.  This 
gives  an  electronic  statistical  weight  of  2  to  the  2E+  state  and  4  to  the  2II  state.16 

The  line  strength  is  the  product  of  the  band  strength  and  the  Honl-London  factor  for  the  tran¬ 
sition  divided  by  the  sum  of  the  Honl-London  factors  for  all  transitions  in  the  band.16 

S  =  (Honl-London)(Band  strength)/(2-$0^)(2S  +  l)(2J+ 1)  (5-16) 

The  Honl-London  factors  have  been  tabulated  for  the  2E-2JI  case.15  Absolute  values  of  the  band 
strength  for  the  y(0,0),  y(l,0),  y(2,0),  and  y(3,0)  bands  have  been  measured.19  For  other  bands  the 
band  strength  may  be  approximated  as  the  product  of  a  Franck-Condon  factor  and  the  square  of  the 
electronic  transition  moment. 

Band  strength  =  qvVRe2  (5-17) 

The  electronic  transition  moments  have  been  reported  as  a  function  of  the  r-centroid  for  the  v’  =  0. 
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I,  2,  and  3  progression.17  These  are  relative  transition  moments  for  each  progression  expressed  is 
R0=(oonstant)xf(r-centroid).  Comparison  of  the  qv.v„R,2  product  and  the  absolute  band  strengths18  al¬ 
lows  the  constant  to  be  evaluated  and  all  band  strengths  to  be  on  an  absolute  scale. 

The  rate  constant  kioi)  is  derived  from  the  ionization  cross  sections  reported  for  v’  =  0.  1.  and 
2. 22  The  Einstein  B  coefficient  is  related  to  the  cross  section  through  the  relation23 

Bi2  =  J^dv/h-V  (5-18) 

<rv  is  the  experimental  ionization  cross  section  at  the  frequency  vQ.  The  integral  is  over  the  laser 
bandwidth  which  is  approximated  as  a  gaussian.  The  resulting  B  coefficient  is  changed  to  units  ap¬ 
propriate  to  energy  density  by  multiplication  by  the  speed  of  light.  NO  has  a  very  high  ionization 
cross  section.  The  rate  constant  for  ionization  may  be  of  the  same  order  of  magnitude  as  the  in¬ 
duced  absorption  and  spontaneous  emission  rate  constants  for  weaker  transitions. 

3.  RESULTS  AND  DISCUSSION 

A.  COMPARISON  WITH  EXPERIMENT 

With  die  connection  made  between  tiie  rate  constants  and  the  spectroscopic  parameters,  the  ion 

yield  for  given  laser  intensity  and  pulse  duration  may  be  calculated.  Figure  5-2  shows  a  plot  of  level 

population  vs.  time  for  the  y-(O.O)  RjjflO.S)  transition  of  NO  with  a  laser  intensity  of  150  mW/cm2. 
bandwidth  of  Icin'1.  The  populations  of  X  and  A  reach  a  steady  state  value  determined  by  the  ratio 
of  statistical  weights  (electronic  and  rotational)  for  the  two  levels.  As  time  increases  the  population  of 
these  two  levels  decreases  due  to  production  of  ions  and  the  spontaneous  emission  from  A  to  other  v" 
and  J"  in  X.  The  population  of  C  and  hence  the  ion  signal  is  not  separately  a  function  of  laser  in¬ 
tensity  and  pulse  duration  but  depends  on  the  fluence,  IAt.  For  the  rest  of  this  paper,  ion  signals 
will  be  reported  as  functions  of  I6t  in  mJ/cm2. 

Zare  et.  a/.14  reported  isopower  spectra  for  a  portion  of  the  y(0.0)  band  at  differing  integrated 
laser  powers  (Fig.  5-3a).  Fig  5-3b  shows  tiie  same  spectra  as  produced  using  the  kinetic  model 
described  above.  For  converting  from  the  reported  integrated  laser  power  to  energy  density,  a 
bandwidth  of  0.7  cm  1  has  been  assumed,  estimated  from  the  resolution  of  rotational  lines.  The  tran¬ 
sition  frequencies  for  the  Q,  -t-  P2,  and  branches  as  determined  by  the  molecular  constants17  do 
not  match  exactly  with  the  positions  in  the  experimental  spectra.  In  the  synthesized  spectra  the 
Qi+P2]  branch  was  shifted  to  give  tiie  same  overlap  as  in  the  experimental  spectra.  The  relative 
agreement  between  the  features  of  the  experimental  and  synthesized  spectra  for  different  integrated 
laser  powers  illustrates  the  effectiveness  of  die  the  kinetic  model. 

The  quantitative  capabilities  of  the  kinetic  model  are  better  exemplified  by  comparison  with  a 
plot  of  ion  signal  vs.  integrated  laser  intensity.  Fig.  5-4a  is  such  a  plot  for  the  R,,(I0.5)  transition 
in  the  y(0.0)  band  as  obtained  experimentally.14  Fig  5-4b  shows  the  ion  signal  dependence  on  laser 
power  for  the  same  transition  as  obtained  from  the  kinetic  model.  Since  in  die  experimental  spectra 
the  R21(10.5)  line  is  not  distinguished  from  the  Q,+P2I(21.5)  the  ion  signal  is  taken  as  the  sum  of 
the  two  transitions. 
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Since  (l  +  l)  MPI  requires  the  absorption  of  two  photons  for  an  ionization  event,  the  dependence 
of  the  ion  signal  on  the  laser  power  should  be  quadratic.  When  the  resonant  ground  to  intermediate 
state  step  is  saturated,  the  ion  signal  depends  only  on  the  rate  of  ionization  and  should  be  linear  with 
respect  to  the  laser  power.  Fig  5-4  shows  that  the  ion  signal  is  neither  quadratic  nor  linear  with 
respect  to  the  laser  power.  To  understand  more  clearly  this  intermediate  behavior  of  the  ion  signal, 
it  is  useful  to  prepare  a  plot  of  the  ln(ion  signal)  vs.  in(integratcd  laser  power).  Fig.  5-5.  Since  the 
ion  signal  is  proportional  to  the  laser  power  to  some  exponent,  S^r.  the  slope  of  a  log  plot  indicates 
the  laser  power  dependence.  The  figure  shows  that  over  a  great  range  the  dependence  is  quadratic, 
and  at  high  powers  goes  to  zero.  The  zero  dependence  on  the  integrated  laser  intensity  means  the 
energy  density  is  great  enough  that  within  the  duration  of  the  laser  pulse  all  molecules  within  the 
beam  are  converted  to  ions.  Clearly  this  would  be  the  ideal  region  to  work  in,  but  it  is  unfor¬ 
tunately  a  few'  orders  of  magnitude  away  from  current  laser  powers.  The  transition  region  between 
the  quadratic  and  zero  dependences  is  where  current  laser  powers  allow  experiments  to  be  performed 
Unfortunately  this  is  also  the  region  where  the  analysis  of  the  ion  signal  is  most  difficult. 

B.  VIBRATIONAL  SUBBAND  RELATIONS 

An  expression  for  the  ion  signal  may  be  written  as  in  equation  5-19  where  S  is  the  ion  signal, 
a  is  an  instrument  constant,  Njo|)  is  the  fraction  of  the  NO  number  density  within  the  laser  beam 
converted  to  ions  in  one  pulse,  and  [NO]  is  the  NO  number  density. 

S  =  aN.JNO]  (5-19) 

The  quantity  of  interest  here  is  the  fraction  of  ions  produced.  Ni0[).  It  has  been  assumed  in 
(l  +  l)  MPI  studies  of  the  vibrational  population  distribution  of  ground  state  NO  molecules6'12  that  at 
saturation  conditions  the  differences  in  bandstrengths  will  no  longer  have  an  effect  on  the  spectra. 
The  form  of  the  spectra  is  determined  only  by  the  differences  in  the  statistical  weights  of  the  X  and 
A  levels  and  the  Boltzmann  distribution  in  the  X  state.  The  ion  signals  should  be  linear  witli  respect 
to  laser  fluence.  and  recording  signals  in  A/B  mode  should  correct  for  fluctuations  in  laser  power. 
This  picture  assumes  a  very  rapid  equilibration  between  the  X  and  A  levels,  an  ionization  rate  much 
less  than  the  induced  absorption  or  stimulated  emission  rates,  and  an  ionization  cross  section  inde¬ 
pendent  of  wavelength  and  intermediate  state  quantum  numbers.  These  conditions  are  not  met  exactly, 
but  with  the  kinetic  model  described  above,  it  is  straightforward  to  calculate  how  the  ion  signal  varies 
for  different  vibrational  bands.  Figure  5-6  shows  Nion  divided  by  the  integrated  pulse  energy  as  a 
function  of  integrated  pulse  energy  for  the  first  four  vibrational  levels  of  NO  at  the  R^OO-.S)  line. 
At  sufficiently  high  laser  fluence.  the  flat  slope  expected  for  linear  dependence  is  observed,  but  the 
fraction  of  ions  produced  in  the  different  vibrational  bands  is  different  at  all  integrated  laser  inten¬ 
sities.  This  means  the  same  ion  signal  measured  at  the  same  laser  power  for  different  vibrational 
bands  corresponds  to  different  ground  state  number  densities.  Clearly  any  measurements  of  the  vibra¬ 
tional  population  distribution  will  need  to  correct  for  this  variation. 

For  measuring  vibrational  populations  using  (1  +  1)  MPI,  we  make  the  following  suggestions. 
Ion  signals  should  be  recorded  in  the  method  outlined  by  Jacobs  and  Zare13.  From  the  curves  of 
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ion  signal  vs.  integrated  laser  intensity,  a  common,  central  integrated  laser  intensity  is  chosen.  Thr> 
bandwidth  of  the  laser  and  the  average  pulse  duration  should  be  well-known  experimentally.  The  frac¬ 
tion  of  the  ground  state  number  density  converted  into  ions  at  these  conditions  is  calculated  for  each 
of  the  vibrational  bands  measured  using  a  model  such  as  described  above  or  perhaps  a  more  compli¬ 
cated  model  if  the  experimental  conditions  necessitate  it.6  A  calibration  curve  of  the  ground  vibrational 
state  is  made  using  known  number  densities.  The  ratio  of  the  fraction  of  ions  produced  in  the 
ground  vibrational  level  to  all  of  the  other  vibrational  levels  is  used  to  correct  the  ion  signals  to  the 
calibration  curve  for  the  ground  state.  For  example,  if  20%  of  the  Y-(0.0)  band  is  converted  to  ions 
but  only  10%  of  the  y*(1.0  band  is,  the  ion  signal  for  the  Y-(l.l)  band  must  be  multiplied  by  two 
to  represent  the  same  percent  of  the  number  density  at  which  the  calibration  curve  was  made. 

C.  CORRECTION  OF  PREVIOUS  KINETIC  RESULTS 

In  reference  II  the  product  vibrational  distribution  of  the  N(4S)+O2-»NO(v=0-7)  was  measured 
using  (l  +  l)  MPI.  The  variations  in  ion  signal  for  the  different  vibrational  levels  was  not  taken  into 
consideration  when  this  work  was  done.  We  will  correct  these  kinetic  results  using  the  methods  out¬ 
lined  here.  The  kinetic  runs  were  recorded  in  A/B  mode,  so  from  a  plot  such  as  was  shown  in 
Fig.  5-6.  we  will  be  able  to  correct  for  the  different  ion  fractions  for  each  of  the  subbands.  Unfor¬ 
tunately  the  data  was  recorded  at  the  P,  bandhead.  While  this  is  the  most  intense  spectral  feature,  it 
is  composed  of  many  transitions.  This  makes  estimating  the  fraction  of  ions  produced  more  difficult 
because  many  lines  must  be  summed.  Figure  5-7a  to  5-7g  is  a  plot  similar  to  figure  5-6  except 
these  are  the  ion  fractions  for  the  P,  bandhead. 

4.  CONCLUSIONS 

A  method  of  relating  ion  signals  measured  in  (l  +  l)  MPI  to  ground  state  number  densities  for 
different  vibrational  bands  has  been  presented. 

The  kinetic  model  outlined  in  this  paper  gives  good  agreement  with  the  experimental  ion  signal 
dependence  on  laser  fluence.  The  good  agreement  with  experiment  suggests  that  the  model  may  be 
used  as  a  method  of  obtaining  absolute  ionization  cross  sections.  Experimental  ion  signal  versus  in¬ 
tegrated  laser  power  curves  may  be  fitted  by  the  model  using  the  ionization  cross  section  as  the  only 
adjustable  parameter. 

Using  normalization,  Franck-Condon  and  our  knowledge  of  the  usual  form  of  linear  surprisals 
along  with  the  approach  outlined  here,  we  should  be  able  to  correct  the  results  for  the 
N (4S) + 02+N0(v = 0-7) + 0(3P)  reaction . 
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FIG.  1.  A)  Schematic  drawing  representing  the  processes  of  induced  absorption,  stimulated  emis¬ 
sion.  spontaneous  emission,  and  ionization  for  a  noncoincident  transition,  a  3-ievel  system.  B)Simi)ar 
schematic  for  a  coincident  transition,  a  4-ievel  system. 
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6.  Stimulated-Emission-Pumping  of  Nitric  Oxide 
6.1.  Introduction 

There  is  much  interest  in  dynamics  of  high  vibrational  levels  of  nitric  oxide.  NO(v=l  to  18) 
is  formed  chemically  in  the  upper  atmosphere,  and  contributes  significantly  to  atmospheric  infrared 
luminosity1.  By  dissociating  NOC1  in  a  pulsed  electrical  discharge,  Deutsch  observed  NO  infrared 
laser  action  in  the  5.84  to  6.43  pm  range,  corresponding  to  the  v=6-ll  fundamental  emission2. 
Experiments  on  an  ultraviolet  NO  laser,  pumping  the  y(0,0)  band  at  227  nm,  have  produced  10  pJ 
pulses  at  237  and  248  nm,  corresponding  to  the  y(0,l)  and  7(0,2)  bands3.  In  order  to  obtain  a 
quantitative  description  of  these  phenomena,  it  is  necessary  to  determine  the  state-specific 
vibrational-to-vibrational  (V-V)  and  vibrational-to-translational  (V-T)  energy  transfer  rates. 
Studies  indicate  that  deactivation  of  NO  (v"=l)  by  a  variety  of  collision  partners  is  extremely 
efficient,  compared  to  species  of  similar  size  and  vibrational  frequency  such  as  CO,  N2,  and 
02.1>4'9  Room  temperature  NO  self-collisions  are  roughly  105  times  more  efficient  in  v=l-»0 
deactivation  than  self-collisional  deactivation  in  these  other  gases,  and  this  efficiency  increases  as 
the  temperature  decreases  below  200K.10  Two  suggestions  have  been  made  to  explain  this 
unusual  behavior.  The  first  of  these,  proposed  by  Nikitin,11  is  that  of  the  two  interaction 
potentials  between  two  v=0  free  radicals,  such  as  NO  with  its  odd  electron,  one  (£g)  is  bound, 
while  the  other  (Iu)  is  antibonding.  The  v=l  +  v=0  surfaces  are  the  same,  except  shifted  up  in 
energy  by  one  vibrational  quantum.  Therefore,  the  binding  surface  of  the  v=l  +  v=0  interaction 
crosses  the  antibonding  surface  of  the  v=0  +  v=0  interaction.  Thus,  vibrational  relaxation  may  be 
viewed  as  collision-induced  intersystem  crossing.  While  this  may  explain  the  highly  efficient 
collisional  deactivation,  further  examination  of  this  scheme12  indicates  that  the  rates  should 
decrease  with  temperature,  contrary  to  observation.  The  second  suggestion13  arises  from  the 
observation  that  the  bonding  well  in  these  interactions  is  approximately  1.6  kcal  mole-1  (560  cm-1, 
which  is  about  2.5kT  at  300K).  In  such  a  case,  a  collision  complex  may  form. 

In  addition  to  self-collisional  deactivation  rates,  the  deactivation  rates  with  O2  are  rather  high 
as  well.  Studies1  have  shown  the  NO(v=l-7)  +O2  collisions  are  10  to  100  times  as  efficient  as 
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NO-noble  gas  collisions  in  inducing  vibrational  relaxation.  These  rates  scale  roughly  with  v2. 
This  is  in  conflict  with  theories  on  V-T  transfer,  such  as  SSH  (Schwartz-Slawsky-Herzfeld), 
which  predicts  the  rates  to  scale  linearly  with  v  14. 

Much  work  has  been  done  to  measure  these  rates.  Most  of  this  work  entailed  collisionally 
exciting  NO,  and  monitoring  the  decay  of  the  infrared  emission.  Methods  which  have  been  used  to 
collisionally  activate  NO  include:  high  energy  electron  bombardment,1*2  electronic-to  vibrational 
energy  transfer  from  mercury  atoms,9  and  gas  discharge15.  These  methods  populate  a  range  of 
vibrational  levels,  rather  than  individual  levels.  Therefore,  extensive  deconvolution  of  the  data  is 
required  in  order  to  retrieve  the  relation  rates  of  individual  quantum  states.  Additionally,  it  is  often 
impossible  to  distinguish  V-V  from  V-T  energy  transfer  rates  in  these  experiments,  leading  to  large 
uncertainties. 

A  number  of  experimental  programs  under  way  in  the  Optical  Physics  and  Infrared 
Technology  Divisions  of  the  Air  Force  Geophysics  Laboratory  require  detailed  state-to-state 
vibrational  relaxation  information  for  NO,  including  the  COCHISE16'19,  LABCEDE1, 
FAKIR17*19,  and  FACELIF20  experiments.  The  observation  of  NO  vibrational  distributions  is 
common  to  all  of  these  experiments.  However,  in  these  experiments,  some  vibrational  relaxation 
occurs  between  formation  and  detection  of  these  distributions.  While  vibrational  deactivation  rates 
are  determined  in  the  LABCEDE  work,  vibrational  excitation  is  not  state-specific,  and  the 
assumptions  which  must  be  made  in  the  deconvolution  result  in  significant  uncertainty  (20-50%  for 
v=2-6).  More  accurate  rates,  especially  for  deactivation  by  O2,  would  improve  the  accuracy  of  all 
of  these  experimental  results.  This  improvement  may  be  realized  by  employing  an  experimental 
method  in  which  one  populates  a  single  vibrational  level  and  monitors  the  decay  of  this  population. 
In  this,  there  is  no  repopulation  from  a  higher  vibrational  level  (cascading  mechanism),  and  the 
analysis  is  considerably  less  complex. 

As  an  example  of  this  need  for  improved  measurement  of  NO(v)  deactivation,  consider  the 
results  from  FACELIF.20  Plots  of  [NO  (v)]  versus  [O2]  show  downward  curvature  at  high  [O2]. 
The  NO  concentration  is  low  enough  in  these  experiments  that  self-collisionai  relaxation  is  not 
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thought  to  introduce  significant  error.  However,  the  O2  reactant  is  in  high  concentration,  so  that 
errors  in  the  rates  of  02-induced  relaxation  processes  can  lead  to  considerable  uncertainty  in  the 
nascent  vibrational  distribution,  and  thus  in  the  state-specific  formation  rates.  This  curvature  has 
been  analyzed  to  estimate  the  rate  for  deactivation  by  O2.  Agreement  between  these  results  and  the 
LABCEDE  results  was  within  experimental  error  (about  30%). 

In  addition  to  bettering  the  accuracy  of  the  i~tes  of  these  important  atmospheric  reactions, 
selective  excitation  of  individual  NO(v")  states  would  allow  direct  investigation  of  interesting 
mechanisms  which  seem  to  be  occurring  in  these  reactions.  Most  notably,  the  formation  of  high-J 
band  heads  in  COCHISE  has  been  attributed  to  the  N(2P)  reaction.  However,  one  may  also 
suggest  that  NO(v)  +  O2  vibrational  deactivation  may  proceed  preferentially  by  vibration-to- 
rotation  energy  transfer  to  high-J  levels.  This  suggestion  is  further  strengthened  by  noting  that  one 
possible  explanation  for  the  extremely  fast  deactivation  is  the  formation  of  an  orbiting  collision 
complex. 

Two  methods  of  selectively  populating  vibrational  levels  can  be  considered.  The  first,  direct 
overton  pumping,  has  been  used  to  study  v=2  relaxation  rates,  including  rotational  and  spin 
relaxation.21  However,  because  overtone  absorption  cross- sections  tend  to  decrease  geometrically 
with  the  change  in  v,  higher  vibrational  levels  cannot  be  populated  directly.  Stimulated  emission 
pumping  (SEP)  was  suggested  as  a  more  effective  method  of  populating  levels  v=6  through  v=16. 

Kittrell  and  coworkers  have  demonstrated  stimulated  emission  pumping  of  iodine.22  The 
"folded  double  resonance"  or  "pump/dump"  technique  utilizes  two  fully  allowed  transitions: 
excitation  to  an  upper  electronic  state  ("pump"),  followed  by  stimulated  emission  to  a  selected 
vibrational  level  in  the  ground  electronic  state  ("dump").  This  technique  is  limited  only  by  Frank- 
Condon  factors,  and  has  been  applied  to  a  number  of  chemical  species,  including  HCN,23 
formaldehyde24-25,  and  p-difluorobenzene.26'29  One  may  tune  the  dump  laser  while  monitoring 
the  fluorescence  from  the  initially  pumped  state.  When  the  dump  laser  wavelength  corresponds  to 
a  resonance,  one  will  see  a  decrease  in  the  fluorescence  due  to  the  depletion  of  the  population  in  the 
initially  excited  state.  This  is  called,  appropriately  enough,  "Fluorescence  Dip  Spectroscopy." 
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This  is  used  more  as  a  diagnostic  for  correctly  setting  the  dump  laser  wavelength,  but  can  be  used 
to  elucidate  the  spectroscopy  of  molecules  which  had  not  been  determined  by  conventional  means 
(e.g.  dispersed  fluorescence  spectroscopy). 

Once  the  population  in  the  target  vibrational  state  has  been  generated,  the  decay  of  this 
population  may  be  followed  using  any  number  of  probe  techniques.  These  include  Laser  Induced 
Fluorescence,  Multiple  Photon  Ionization,  and  Infrared  Fluorescence  or  Absorption.  If  LIF  is 
used,  the  technique  is  commonly  called  "Triple  Resonance",  in  which  the  target  level  population  is 
mapped  out  as  the  probe  laser  is  delayed  in  time. 


6.2  Evaluation  of  Available  Pump  Sources 
6.2.1.  SEP  of  NO  using  the  ArF  excimer  pump 

Our  initial  experiments  used  an  ArF  excimer  laser  as  the  pump  source.  NO  absorbs  193  nm 
light  via  high  J  transitions  of  the  B2I"I(v'=7)  X2n(v'=0)  [0(7,0)]  band30.  Because  of  the 

broad-band  nature  of  the  ArF  excimer  (approximately  100  cm-1),  many  rotational  states  are 
populated  via  the  transitions  P(25-30)  and  R(28-33).  One  needs  to  stimulate  emission  from  all  of 
these  levels  in  each  laser  shot  for  two  reasons.  First,  if  only  one  rotational  level  is  depopulated, 
only  a  slight  decrease  in  the  fluorescence  may  be  observed,  rendering  it  very  difficult  to  locate  the 
fluorescence  dip  signal.  Second,  those  molecules  remaining  in  the  upper  state  are  then  free  to 
spontaneously  radiate,  thus  jeopardizing  the  selectivity  of  the  process.  This  is  particularly 
hazardous  in  the  NO  B<-X  transitions,  where  the  Franck-Condon  distribution  is  quite  broad. 
Therefore,  in  order  to  depopulate  all  of  these  levels,  it  is  necessary  to  use  a  broad  band  source  as 
the  dump  laser.  As  the  output  ot  tiiis  laser  could  not  be  frequency  doubled,  we  were  constrained  to 
visible  transitions,  with  target  levels  from  v"=19  to  23. 


The  specifications  of  the  Lumonics  Excimer  Laser  limited  the  effectiveness  of  this  approach. 


Because  the  ArF  excimer  is  broad-band,  the  spectral  brightness  is  low.  The  output  beam  is  also 
quite  large  and  divergent,  necessitat*  .  the  use  of  apertures  and  spatial  filters,  which  lead  to 
significant  loss  of  power,  Finally,  this  f  ed  frequency  laser  cannot  be  tuned  to  the  band  head,  and 
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therefore  accesses  only  high-J"  levels  at  the  tails  of  the  thermal  rotational  population  distribution. 
We  had  intended  to  use  isotopically  substituted  15N180  to  shift  the  maximum  of  the  rotational 
distribution  into  better  resonance  with  the  laser.  This  technique  has  been  used  successfully  in  our 
laboratory,  in  which  isotopically  substituted  methane(13CD4)  is  used  to  bring  a  selected 
rovibrational  transition  into  resonance  with  a  fixed-frequency  CO2  laser  line.31*32  In  this  case,  as  is 
shown  below,  such  isotopic  shifting  was  not  sufficiently  effective  to  overcome  the  other  inherent 
problems. 

In  order  to  quantify  the  limitations  of  the  ArF  excimer  approach  to  the  stimulated  emission 
pumping  of  NO,  we  can  compare  it  to  the  work  done  in  Professor  A.E.W.  Knight's  laboratory  at 
Griffith  University  (Queensland,  Australia)  on  I2  in  a  free  jet  expansion.  This  way,  we  can  obtain 
a  good  estimate  of  the  signal-to-noise  ration  which  may  be  expected,  based  on  experience. 

In  these  experiments,  the  output  from  one  Nd:YAG  laser  (Molectron)  was  split,  and  used  to 
pump  the  two  dye  lasers  (Molectron  and  Lumonics)  which  were  used  in  the  SEP  pump-dump  step. 
Another  Nd:YAG-pumped  dye  laser  (Molectron)  was  used  as  the  probe.  Each  of  these  provided 
10  ns  pulses  with  a  few  mJ  in  the  visible.  The  (20,0)  transition  near  560  nm  was  pumped,  and 
fluorescence  from  the  v'=20  level  was  dispersed.  We  set  the  monochromator  on  the  (20,4)  dump 
transition  at  586  nm,  then  tuned  the  dump  laser  until  scattered  light  at  this  wavelength  was 
maximized.  The  PMT  was  equipped  with  a  band-pass  filter  (^=532  nm,  Ak=10  nm  FWHM  to 
protect  it  from  scattered  light,  and  also  from  fluorescence  due  to  pump  only.  With  this  mounted  on 
the  vacuum  chamber,  we  tuned  the  probe  laser  to  the  (24,4)  transition  at  555  nm,  and  observed  the 
(24,0)  fluorescence  at  532  nm.  The  triple-resonar.ce  signal  was  obtained  quite  easily  (S/N  >  25). 
However,  because  lifetimes  in  the  X3Z  state  are  much  longer  than  in  the  B3I1  state,  molecular  flow 
out  of  the  detection  region  dominated  any  collision  dynamical  effects. 

First,  consider  the  pumping  rate  for  the  (20,0)  transition  at  560  nm.  The  fluorescence 
lifetime  is  0.66  (ts.33  Using  the  relations  A=r]  and 


we  obtain  a  value  for  the  Einstein  B  coefficient  for  I2  B3I1  «-  X3Z  vibronic  transitions,  for  light 
intensity  in  erg  s*1  cm*2: 

B1  =  5.3  x  107sec  g'1 

Typically,  with  the  second  harmonic  of  the  Nd:YAG  laser  pumping  two  dye  lasers,  10  ns 
pulses  of  1  mJ/pulse  and  Av=l  cm-1  are  obtained.  The  pump  rate  ki2PumP  =  BI/Av.  For  these 
parameters  ki2PumP  =  7.1  x  1011  s_1,  five  orders  of  magnitude  faster  than  the  spontaneous 
fluorescence  rate.  Therefore,  the  pump  transition  is  saturated,  and  50%  of  the  population  is 
moved. 

The  stimulated  emission  rate,  or  ki2^umP  is  obtained  by  correcting  for  Franck-Condon 
factors34  and  transition  frequencies.  For  dumping  into  v"=4  at  586  nm,  the  rate  is  9  x  1011  s'1. 
Therefore,  this  transition  is  also  saturated.  Thus,  25%  of  the  initial  population  in  the  ground  state 
is  moved  into  the  selected  target  vibrational  level. 

Only  one  rotational  level  may  be  accessed  with  this  narrow-band  laser.  The  laser  is  tuned  to 
transitions  from  J"max-  At  X/D=5,  T=12K,  and  the  rotational  degrees  of  freedom  have  not  yet 

frozen  out.  The  rotational  constant  of  I2  is  .04  cm’1-35,  resulting  in  Jmax  =  10.  0roi  =  .054K,  so 
qrot  =  222.36  The  fraction  of  molecules  in  J"=10  which  are  available  for  pumping,  Pjmax,  is  0.06. 

The  actual  number  density  of  I2  molecules  in  the  jet,  under  typical  operating  conditions  (2% 
12  in  Ar,  p0  =  1  atm,  X/D=5),  Pj2>  is  1017  molec  cm-3.  Using 

pl2,v"=4  =  pl2  X  PJmax  X  P(3umP  X  PPumP 
we  find  the  population  of  molecules  in  the  target  vibrational  levels  to  be: 

pl2  V,I=4  =  3  x  1013  molec  cm-3 

The  corresponding  parameters  for  NO  P(7,0)  are  A  =  2.5  ps,37  Bi=5.3  x  105  sec  g'1,  and 
the  laser  pulses  are  10  ns  long,  with  about  100  pJ/pulse,  Av  =  100  cm-1.  The  pulse  energy  is 
lower  than  one  would  expect  because  of  the  difficulties  of  shaping  and  steering  the  excimer  beam. 
Following  an  analogous  calculation  to  that  above  for  I2,  we  determine  that  the  pumping  rate  for  the 
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NO  P(7,0)  transition  is  1.7  x  106  s'1.  This  must  be  compared  to  the  inverse  of  the  zero  pressure 
lifetime  (including  spontaneous  fluorescence  and  predissociation)  of  0.331  |j.s,30  giving  a  rate  of  3 
x  106  s'1.  Since  these  rates  are  of  the  same  order  of  magnitude,  we  must  delve  further  to  determine 
the  actual  fraction  of  molecules  pumped  under  these  conditions. 

Toward  this  end,  we  can  consider  this  as  a  three- level  system,  as  shown  in  Figure  6-1.  The 
rate  equations  for  this  system  have  been  solved.38  Including  the  simplification  that  ki  =  k.i  [Boi= 
Bio],  the  number  densities  in  all  three  levels  are  described  by  the  following  equations: 

N0  =  ~-2— kl  e~mit  +  — l-~ml  e~m2t 

m2  “  mi  m2  -  mi  (6-3a) 

Nj  = — ^ — fe-mll_e-m2t] 

m2  ~mi  1  J  (6-3b) 


Nout  =  H—.ml  e~rn2t - E2_e-mit 

m2  “mj  m2-m1  (6-3c) 


where: 

mi  =  l/2[(2ki  +  k2)  -  (4ki2  +  k22) 1/2] 
m2  =  l/2[(2ki  +  k2)  +  (4ki2  +  k22)1/2J 
m2  -  mj  =  (4kj2  +  k^)1/2 

Figure  6-2  shows  the  results  of  these  calculations.The  lower  plot  is  for  the  more  general 
case,  indicating  that  the  pumping  rate  must  be  approximately  two  orders  of  magnitude  greater  than 
the  rate  of  removal  from  the  Rabi  cycle  in  order  to  ensure  saturation.  The  upper  plot  is  for  this 
case,  with  parameters  particular  to  the  Lumonics  Excimer  laser  (1  mm2  spot  size,  Av=100  cm-1, 
and  the  laser  pulse  is  approximated  by  a  10  ns  square  pulse).  From  this  we  see  that,  with  100 
joJ/pulse,  only  about  1%  of  the  population  is  moved  into  B2fl  (v-7). 
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Now  consider  the  "dump"  transition.  kdump  may  be  obtained  by  correcting  kpUmp  for 
Franck-Condon  factors,39  transition  frequency,  and  laser  power,  in  the  same  way  as  was  done  for 
12  above.  The  BBDL  delivers  about  300  (iJ/pulse,  which  gives  a  dump  rate  for  the  p(7,23) 
transition  at  610  nm  as  kdump  =  -84  x  106  s*1.  Approximately 
f  dt  ~  At  and  Jb23  (oo-co0)  I(o(G>-tt0)dco  -  kdump 
and  using  the  expression  from  Kittrel,  et  al.:22 

n2(0Vt)  =  (l/2)n2(03,0)  x  [1  +  exp(-2kdumpAt)J  (6-4) 

then:  n2(co=  10ns)  =  G.92  n2(co,'t) 

and  only  8%  of  the  population  is  moved  from  B(v’=7)  into  the  target  vibrational  state  in  X. 

We  follow  a  similar  calculation  to  that  above  to  determine  the  populations  in  the  relevant 
rotational  states.  The  spectroscopic  constants  are  found  in  Huber  and  Herzberg.35  If  the  14N160 
is  used,  J=25-28  may  be  pumped.  The  total  fraction  in  these  states  is  about  .012.  These  high 
rotational  levels  are  pre-dissociative,  giving  rise  to  higher  oxides  of  nitrogen  in  recombination 
reactions.  In  order  to  freeze  out  these  impurities,  the  sample  must  be  held  over  a  liquid  nitrogen 
trap.  NO  has  a  vapor  pressure  of  about  100  mTorr  at  77K.40  At  100  MTorr,  pN0  =  3.5  x  1015 

molec  cm-3,  and 

14N160:  pv"=23  =  3.5  x  1010  molec  cnr3 
The  isotope  15NlgO  helps  somewhat,  giving  access  to  J"  =  15-18  ,  so 
15N180:  pv"=23  =  2.1  x  1011  molec  cm-3 

The  probe  and  detection  transitions  for  both  experiments  are  similar  in  that  they  have  large  Franck- 
Condon  factors,  the  lasers  used  are  of  similar  bandwidth  and  power,  and  that  a  band-pass  filter  is 
used  to  separate  the  desired  triple  resonance  fluorescence  from  the  fluorescence  due  to  the  pump 
only  step.  However,  we  can  see  that  the  number  density  in  the  14N160  target  level,  even  in  the 
best  case,  is  a  factor  of  103  lower  than  that  for  I2  SEP  work.  Empirically,  we  obtained  a  signal-to- 
noise  ratio  in  the  I2  work  of  about  25.  Scaling  this  to  the  NO  work,  we  would  expect  a  signal-to- 
noise  ratio  of  .025.  Even  with  the  heavy  isotope,  S/N  is  still  below  1 ,  rendering  it  nearly 
impossible  to  find  the  triple  resonance  signal  at  all 
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6.2.2  SEP  of  NO  using  a  Two-Photon  Pump 

Given  the  limitations  of  using  the  ArF  excimer  laser  as  a  pump  source,  we  must  consider 
available  tunable,  narrow-band  sources  for  exciting  NO.  The  7(0,0)  transition  at  226  nm  is  well 
known,  and  has  been  the  basis  of  a  number  of  experiments.  However,  the  y-bands  follow  a  Av=l 
propensity  rule,  due  to  the  similarity  between  the  X(2I7)  and  A(2E)  potentials.  As  a  result, 
A(2£)(v'=0)  has  small  Franck-Condon  overlap  with  any  high  vibrational  levels  of  X(2n),  limiting 
our  ability  to  populate  any  of  these  levels  of  interest.  Higher  vibrational  levels  of  the  available 
excited  electronic  states  (A,B,  and  C)  have  more  advantageous  Franck-Condon  overlaps,  but 
transitions  into  these  states  are  205  nm  and  below,  beyond  the  range  of  available  doubling  crystals. 
Tunable,  narrow-band  VUV  radiation  has  been  generated  by  third-harmonic  generation  in  a 
strontium-xenon  mixtures,41  by  four- wave  mixing  in  mercury  vapor,42  or  by  frequency-doubling 
and  Raman  shifting  the  output  of  a  single  dye  laser.43  The  pulse  energies  obtained  with  these 
systems  are  on  the  order  of  several  (Jj,  which  is  sufficient  for  one-laser  excitation  experiments,  but 
is  not  enough  to  actually  significantly  populate  the  upper  state.  One  may  suggest  an  excimer  laser 
in  which  one  mode  is  selected  by  prisms  in  an  oscillator  cavity,  and  is  amplified  in  another  cavity. 
This  would  give  high-power,  narrow  bandwidth  VUV  radiation,  and  has  been  successfully 
applied  to  KrF  output  near  248  nm.  However,  experience44  has  shown  that  extension  of  this  to 
ArF  output  near  193  nm  is  unreliable  at  best.  Recently,  Sum  Frequency  Generation  in  barium 
borate  (BBO)  nonlinear  crystals  has  been  reported.  In  this,  248.5  nm  (obtained  by  frequency 
doubling  the  output  of  a  dye  laser,  not  from  a  KrF  excimer  laser)  and  800  nm  light  are  mixed  to 
give  10-100  pJ/pulse  near  191  nm  45  This  may  be  promising,  despite  the  limitation  that  the  full 
triple  resonance  scheme  would  then  require  four  dye  lasers,  three  of  which  are  pumped  by  the 
same  source. 

While  generating  high-power  VUV  radiation  remains  elusive,  progress  has  been  made  in 
lasers  operating  in  the  visible  and  near  UV.  Excimer-pumped  dye  lasers,  delivering  pulse  energies 
of  tens  of  ml,  have  introduced  the  possibility  of  nearly  saturating  a  two-photon  transition.  As  will 
be  shown  below,  pumping  transitions  at  the  two- photon  equivalent  of  VUV  frequencies,  followed 
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by  dumping  with  one  UV  photon,  is  an  efficient  method  of  generating  significant  populations  in 
selected  vibrational  levels  of  the  ground  electronic  state  of  NO. 

As  before,  let  us  calculate  the  population  in  the  target  vibrational  level  we  may  obtain  using 
this  method.  Several  values  for  the  two-photon  absorption  cross-section  for  the  y(0,0)  [R22  +  S12 
(J"=9.5)]  transition,  integrated  over  the  Doppler  line- width  of  0.1  cm-1,  have  been  reported  in  the 
literature.46*49  These  values  range  from  3  x  10'48  to  3  x  10'49  cm4  sec.  In  this  order-of- 
magnitude  calculation,  let  us  simply  use  1048  cm4  sec.  So,  for  a  laser  pulse  with  a  band-width 
equal  to  the  Doppler  width,  the  rate  of  pumping  is: 


_  2jioI2 
pump  "  (hl))2 


(6-5) 


If  we  pump  the  y(3,0)  transition  near  391  nm  under  typical  experimental  conditions  (0.05  mm2 
spot  size,  etalon- narrowed  bandwidth  of  Av=0.04  cm4),  we  obtain  a  pumping  rate  per  mJ  of  laser 
pulse  energy: 

kpump/mJ  =  2.32  x  106  mJ'2sec4. 

This  vibrational  level  is  not  above  the  dissociation  limit,  so  the  only  method  of  removal  from 
the  Rabi  cycle  is  via  spontaneous  fluorescence.  The  fluorescence  lifetime  is  200  ns.50'53  which 
gives  kout  equal  to  5  x  106  sec'1.  Using  the  same  method  as  depicted  in  Figure  6-1,  with  the 
exception  that  ki  and  k.i  are  two-photon  processes,  we  may  calculate  the  fraction  in  each  of  the 
levels  after  a  10  ns  pulse,  as  a  function  of  laser  pulse  energy.  The  results  of  this  calculation  are 
given  in  Figure  6-3. 


6.3  Experimental 

6.3.1.  Two-Photon  Pumping 

This  calculation  of  two-photon  pumping  efficiency  was  experimentally  tested.  The  etalon- 
narrowed  output  of  a  XeCl  excimer-pumped  dye  laser  (Lambda  Physik  EMG201/FL3002)  is 
focussed  into  a  cell  containing  100  mTorr  neat  NO.  Fluorescence  is  filtered  from  scattered  laser 
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light  (Corion)  and  imaged  onto  a  photomultiplier  tube.  The  laser  is  pressure  scanned  with 
nitrogen,  and  a  small  portion  of  the  two-photon  laser  excitation  spectrum  is  given  in  Figure  6-4 
with  assignments  after  Zacharias,  et  al..32  The  laser  is  then  set  on  the  Pn(9  1/2)  transition.  A 
back  reflection  impinging  on  a  photodiode  is  used  to  measure  relative  laser  energy,  and  the  beam  is 
dumped  into  a  joule  meter  (Laser  Precision)  for  absolute  laser  energy  measurements.  The  laser  is 
attenuated  by  a  series  of  glass  microscopy  cover  slips  (each  attenuates  by  about  10%),  and  the 
pulse  energy  and  fluorescence  intensity  are  simultaneously  measured  over  a  wide  range  of 
energies.  The  output  of  the  PMT  and  the  photodiode  are  integrated  (EG&G/PAR)  and  collected 
on  a  laboratory  computer  (IBM  PC). 

The  data  are  then  transferred  to  a  mainframe  computer  (VAX)  for  analysis.  Laser  pulse 
energy  is  scaled  to  the  absolute  values  given  by  the  joule  meter.  The  maximum  signal  intensity  if 
then  scaled  to  the  level  given  by  the  calculation.  As  shown  in  Figure  6-5,  there  is  good  agreement 
between  experiment  and  the  functional  form  of  pumping  efficiency  versus  pulse  energy,  even  with 
this  simple  scaling  scheme,  throughout  the  range  of  energies  obtained.  Therefore,  we  may  be 
confident  of  the  veracity  of  the  scheme  depicted  in  Figure  6-3. 

6.3.2.  "Two  minus  One"  SEP 

With  reasonably  efficient  two-photon  pumping  of  the  7(3,0)  Pn(9  1/2)  transition,  we  may 
now  include  the  second  step,  the  dump  of  population  into  the  target  vibrational  level  of  X(2n) 
state.  For  this,  two  cells  containing  400  mTorr  NO  (1.4  x  1016  molecules  cm'3)  are  used.  The 
pump  beam  is  focussed  into  both  (400  mm  focal  length),  and  the  LIF  signals  from  these  are 
differentially  amplified  and  balanced  to  give  a  zero  baseline.  In  order  to  maximize  the  pumping 
efficiency,  the  pump  beam  was  not  split,  but  rather  focussed  consecutively  into  these  two  cells. 
With  a  typical  pump  energy  of  1.0  mJ,  our  calculation  above  indicates  that  approximately  2%  of 
the  molecules  in  the  focal  region  are  pumped.  The  dump  laser  output  (Lambda  Physik 
EMG200/FL2002,  operating  with  Coumarin  460)  is  frequency  doubled  to  give  200  joJ  pulses  near 
228  nit.,  the  7(3,4)  band.  This  is  optically  delayed  by  10  ns,  then  is  loosely  focussed  (600  nm 
focal  length)  and  directed  colinear  and  counter  propagating  to  the  pump  beam  through  one  of  the 
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cells.  It  is  prevented  from  entering  the  other  cell  with  a  glass  microscope  slide.  The  output  of  the 
differential  amplifier  is  integrated  and  averaged  as  before.  As  the  dump  laser  is  tuned  into 
resonance  with  allowed  rovibronic  transitions  in  the  y(3,4)  band,  we  observe  the  dips  in  the 
differential  fluorescence  signal,  corresponding  to  depletion  of  the  A(2E)(v'=3,  J=  8  1/2)  rovibronic 
level  population.  A  schematic  of  the  relevant  levels  and  transitions  is  given  in  Figure  6-6,  and  the 
experimental  fluorescence  dip  spectrum  in  Figure  6-7.  These  dips  correspond  to  a  depletion  of 
about  20%  of  the  population  in  the  upper  state.  According  to  Equation  6-4  above,  this  transition 
should  be  saturated  under  these  conditions,  which  would  give  50%  depletion.  The  discrepancy  is 
most  probably  because  the  dump  beam  is  focussed  more  tightly  than  is  optimal.  If  the  dump  beam 
does  not  fully  cover  the  volume  defined  by  the  pump  beam,  it  will  not  spatially  access  all  the 
initially  pumped  molecules. 

The  fraction  of  population  in  the  initial  v"=0,  J=  9  1/2  is  7.5%  at  300K.  Thus,  we  find  the 
number  density  of  molecules  in  our  target  vibrational  level  to  be: 

PNO’  v"=4  =  PNO  x  x  ^  pump  x  Pdump 
=  2.2  x  1012  molecules  cm*3 

6.4  Discussion 

These  experiments  indicate  that  2-1  SEP  is  an  effective  way  to  populate  single  vibrational 
levels  in  the  ground  electronic  state  of  NO.  It  is  desirable  to  distribute  the  XeCl  excimer  pump 
beam  between  the  pump  and  dump  lasers,  such  that  the  efficiency  for  each  process  is  optimized. 
This  was  not  possible  in  the  above  experiments,  due  to  laboratory  constraints.  Despite  this,  we 
aheady  achieved  an  order  of  magnitude  improvement  in  our  stimulated  emission  pumping 
efficiency  over  the  ArF-pumping  approach.  Let  us  assume  that  we  have  divided  the  XeCl  excimer 
laser  so  that  it  gives  us  6  mJ/pulse  near  390  nm,  to  be  used  for  the  two-photon  pump.  Our 
calculations  then  indicate  that  the  transition  is  nearly  saturated  and  40%  of  the  molecules  may  be 
pumped.  This  will  not  only  increase  the  number  density  of  the  molecules  pumped,  but  will  also 
improve  the  signal-to-noise  ratio  in  the  fluorescence  dip  spectrum.  Near  saturation,  signal  levels 
are  not  as  sensitive  to  fluctuations  in  laser  power. 
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By  insuring  enough  power  for  the  pump  step,  there  is  not  much  of  the  XeCl  excimer  pump 
beam  left  for  the  dump  laser.  Fortunately,  the  one-photon  y-band  transitions  may  be  saturated 
with  only  a  few  tens  of  pJ  of  pulse  energy.  As  indicated  above,  the  A2I(v)  have  lifetimes  on  the 
order  of  200  ns,  which  gives  rise  to  an  Einstein  B  coefficient  of  5  x  106  s  g*1.  For  50  pJ  pulses, 
focussed  to  0.5.5  mm2  spot,  and  with  a  bandwidth  of  0.2  cm*1,  we  obtain  k  dump  of  7.5  x  109 
sec*1.  We  refer  again  to  the  equation  given  by  Kittrell,  et  al.,22  which  gives  us  a  dumping 
efficiency  of  50%. 

The  initially  pumped  state  in  this  double  resonance  process  is  A2Z  (v'=3,  J -8  1/2).  Unlike 
the  high  rotational  states  of  B2n(v'=7),  dissociation  is  not  a  problem  for  these  low  rotational  levels 
of  A(v'=3).  Therefore,  it  is  not  necessary  to  maintain  the  sample  over  a  liquid  nitrogen  trap  to 
freeze  out  recombination  products,  so  higher  pressures  may  be  used.  At  a  modes  pressure  of  1 
Torr  (3.5  x  1016  molecules  cm*1),  we  obtain  for  the  number  density  of  molecules  in  the  target 
vibrational  level  of  the  ground  electronic  state: 

PNO  v"=4  1014  molecules  cm*3 

Although  this  number  is  accurate  only  to  the  order  of  magnitude,  it  is  obvious  that  the 
molecular  densities  in  the  target  vibrational  level  may  be  obtained  with  this  method  which  are 
significantly  higher  than  those  which  may  be  expected  for  pumping  with  a  193  nm  ArF  excimer 
laser.  Also,  a  wide  variety  of  v"  levels  may  be  populated  by  "2  minus  1"  stimulated  emission 
pumping.  For  example,  A2L  (v'=4)  has  good  Franck-Condon  overlap  with  X2I1  (v"=3,5,6,ll, 
and  12).  Nearly  all  other  v"  levels  of  interest  may  be  accessed  via  v'  =  2  and/or  3.  State-specific 
rates  for  collisional  vibrational  relaxation  of  this  important  atmospheric  molecule  may  now  be 
addressed.  Recently,  stimulated-emission  pumping  of  NO  via  the  B  state  has  been  demonstrated 
using  a  tunable,  narrowband  argon  fluoride  laser54,  in  which  the  laser  output  could  be  tuned 
confollably  across  the  excimer  transition  gain  curve.  This  also  affords  the  possibility  of  collisional 
relaxation  measurements  in  high  vibrational  measurements  of  NO. 
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Figure  6-2  .  193  nm  pumping  efficiency  under  sub-saturation  conditions 
Upper  Trace:  Specific  case,  showing  fraction  in  the  upper  level  after  a  10  ns 
pulse  vs.  pulse  energy,  using  parameters  particular  to  the  Lumonics  Excimer  Laser  (see  text). 

Lower  Trace:  General  case,  showing  fraction  in  each  level  after  a  10  ns  pulse 
vs.  rationa  of  rate  into  over  rate  out  of  the  upper  state,  k,/k0(Jt. 
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Figure  6-3.  Two-photon  pumping  efficiency  under  sub-saturation  conditions 
Fraction  of  molecules  in  each  of  the  levels  shown  in  Figure  1-i 
Dashed  line:  NQ 


Figure  6-4.  NO  r(3,0)  Two-Photon  Excitation  Spectrum 
Assignments  are  after  Zacharias,  et.  al.52 
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Figure  6-5.  Experimental  and  calculated  two-photon  pumping  efficiency 
Dots:  Experimental  points 
Solid  line:  Calculated  efficiency 


-photon  pump  NO  7(3,0)  Pj(9  1/2)  (391.5  nm),  dump  7(3,4) 
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7.1  Recent  Developments  and  Future  Applications 

7.1  REMPI  Detection  of  N2(X1I+g) 

While  making  measurements  of  N(^S)  MPI  signals,  a  number  of  molecular  bands  were 

apparent  in  the  same  spectral  range.  These  bands  were  quenched  only  very  slowly  on  addition  of 
either  02  or  NO.  They  also  became  more  prominent  as  the  mole  fraction  of  molecular  nitrogen 

was  increased.  The  band-head  positions  coincide  with  the  reported  band-head  wavelengths*  for 
the  Lyman-Birge-Hopfield  bands  of  N0  (a*n  ->  X*X  +).  The  single  photon  wavelengths  of 

z  g  g 

the  observed  transition  band-heads,  together  with  their  assignments,  are  listed  in  Table  7.1. 

4 

The  observed  spectra  for  the  (9,17)  and  (13,20)  bands  are  shown  in  Fig.  7.1.  The  N(  S) 
quartet  is  visible  to  the  red  side  of  these  bands.  This  N2  transition  is  permitted  only  by  magnetic 

dipole  and  electric  quadrupole  selection  rules.  This  gives  S  and  O  electric  quadrupole  branches 

(AJ  =  ±2)  and  P,  Q.  and  R  electric  quadrupole  plus  magnetic  dipole  branches  (AJ  =  ±1).  Work  is 

currently  undei  way  to  fit  these  spectra  using  the  kinetic  model  desribed  in  Section  5.  This  will 

enable  quantitative  comparison  of  relative  ground  state  populations  to  be  made  for  high 
vibrational  levels  in  the  ground  state  of  N2.  MPI  spectra  of  N2(X)  have  been  reported 

previously  using  (3+1)^  and  (2+2)^  excitation  schemes.  The  spectra  here  are  obtained  using 
a  (1+1)  excitation  scheme  and  might  therefore  be  expected  to  provide  a  more  sensitive  detection 
method.  A  measurement  of  the  power  dependence  of  the  (9,17)  band-head  was  made  and  the 

results  are  shown  in  Fig.  7.2.  The  power  dependence  of  the  MPI  signal  varied  between  1.6  and 

1.2  over  the  laser  power  range  investigated.  This  suggests  that  the  (a*FI  -»X*E  +)  transition  is 

5  o 

becoming  partially  saturated. 

We  have  been  successful  in  detecting  N2  in  highly  vibrationally  excited  states  (v"  =  15  to 

21)  in  active  nitrogen  flows.  It  is  to  be  expected  that  other  vibrational  levels  will  be  accessible  at 
longer  wavelengths  provided  that  sufficient  energy  is  available  to  reach  the  ionization  limit  using 
only  two  photons,  and  that  the  Franck-Condon  overlap*  is  favorable.  Lower  vibrational  levels 
generally  require  progressively  shorter  laser  wavelengths  and  are  therefore  not  accessible  using 
this  scheme. 
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7.2.  Identification  of  Atomic  and  Molecular  Species  Generated  in  Nitrogen 
Afterglows  with  SF^ 

Attempts  were  made  to  enhance  N-atom  dissociation  in  N(^D)+0  quenching  experiments 
(Section  2)  by  flowing  small  amounts  of  SF^  through  the  microwave  discharge.  While  this 

successfully  increased  both  the  N(^S)  and  N(^D)  concentrations,  it  also  resulted  in  the 
appearance  of  atomic  and  molecular  features  in  the  spectral  region  of  interest.  The  observed 
spectrum  is  shown  in  Fig.  7.3.  The  strong  atomic  line  corresponds  to  the  (1+1)  MPI  transition^ 

S[5p  3Pj]  +-  S[3p4  3P2]  (7-1) 

of  atomic  sulphur.  The  splitting  of  the  upper  J-level  components  is  not  resolved.  The  molecular 
spectrum  corresponds  to  the  (1+1)  MPI  B^E+-+X^Ily2  transition  of  NS^.  The  X^F I3/2 

branch  will  be  apparent  some  220  cm'^  to  the  red  although  this  has  not  been  confirmed 

7 

experimentally.  No  other  band  positions  are  reported  for  this  transition  . 

7.3.  FACELIF  Modifications 

A  number  of  modifications  to  the  FACELIF  apparatus  have  been  implemented  in  the  light 
of  some  of  the  limitations  encountered  in  the  experiments  described  in  the  previous  sections  of 
this  report.  These  improvements,  together  with  their  motivation,  are  discussed  briefly  here.  The 
two  main  aspects  of  these  modifications  are,  firstly,  a  substantially  longer  quartz  flowtube  and, 
secondly,  a  movable  MPI  detector. 

The  30  cm  pyrex  flowtube  has  been  replaced  by  1.2  m  length  high  grade  quartz  tubing  of 
the  same  diameter.  This  has  two  distinct  advantages  over  the  previous  arrangement.  The  first  of 
these  concerns  [N]  measurements  for  faster  flows,  made  using  NO  titrations.  As  desribed  in 
Sec.  2,  in  order  to  obtain  accurate  titration  end-points,  sufficient  flow  time  must  be  allowed  for 
the  N  +  NO  reaction  to  go  to  completion.  The  use  of  this  longer  flowtube  will  enable  these 
titration  measurements  to  be  taken  downstream  of  the  MPI  detector,  thus  extending  the  temporal 
range  of  the  titration  technique. 
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The  longer  flowtube  arrangement  also  provides  ample  room  for  the  installation  of  a 
spectrometer.  A  linear  motion  system  has  been  purchased  which,  when  installed,  will  allow  the 
spectrometer  to  be  accurately  translated  along  the  flowtube  axis.  The  spectrometer  will  provide  a 
useful  additional  diagnostic  tool  for  analysing  the  temporal  evolution  of  the  afterglow. 

The  second  major  modification  is  the  installation  of  a  movable  MPI  detector.  With  the 
previous  set-up,  the  grids  remained  fixed  with  the  amplifier  circuit  attached  directly  to  the  feed¬ 
throughs.  The  new  detector  has  the  amplifier  attached  to  the  grids  inside  the  flowtube.  The 
whole  assembly  may  be  translated  by  means  of  a  1/4"  stainless  steel  tube  which  also  serves  as  a 
guide  for  the  electrical  connections. 

The  detector  incorporates  a  three  grid  system  which  subtracts  the  background  current  in 
the  absence  of  the  laser  pulse,  thus  effectively  eliminating  detection  problems  emanating  from 
ions  present  in  the  flow.  The  laser  beam  is  directed  straight  through  the  wall  of  the  flowtube. 
Despite  the  lack  of  Brewster  angle  windows,  the  laser  power  is  not  significantly  attenuated  by  its 
passage  through  the  flowtube.  MPI  spectra  have  been  successfully  obtained  using  the  new 
FACELIF  configuration. 
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Table  7.1  Observed  and  calculated  bandhead  po.  '->r'  for  the  a*rig->X*X„+  transition  in 


N 


2* 


WaYfileneth/nru 

jtl  ill  Observed..  .Calculated 

7  15  207.71  207.66 

8  16  209.63  209.60 

9  17  211.49  211.45 

10  18  213.37  213.35 

12  19  209.13  209.03 

13  20  211.06  211.00 

14  21  213.01  212.98 
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Figure  7.1.  N-(a  EL— )  spectra  showing  (9,17)  and  (13,20)  bands.  Flow  conditions 

~  5  6 

were  10%  N~  in  He  at  a  total  pressure  of  3.2  Torr. 
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figure  7.2.  Laser  power  dependence  of  N2  (9,17)  band.  The  dependence  ranges  from  1.6  to  1.2 
at  higher  power. 
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Figure  7.3.  MPI  spectrum  of  the  region  around  the  N(^D)  transition  in  the  presence  of 
showing  the  S(3P)  line  and  the  NS(B2X+-»X2ril/2)  band. 


